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Abstract. Given a random quantum state of multiple (distinguishable or indistinguishable) 
particles, we provide an algorithm, rooted in symplectic geometry, to compute the joint proba- 
bility distribution of the eigenvalues of its one-body reduced density matrices, and hence some 
associated physical invariants of the state. 

As a corollary, by taking the support of this probability distribution, which is a convex 
polytope, we recover a complete solution to the one-body quantum marginal problem, i.e., 
the problem of characterizing the one-body reduced density matrices that arise from some 
multi-particle quantum state. In the fermionic instance of the problem, which is known as 
the one-body A''-representability problem, the famous Pauli principle amounts to one linear 
inequality in the description of the convex polytope. 

We obtain the probability distribution by reducing to computing the corresponding distribu- 
tion of diagonal entries (i.e., to the quantitative version of a classical marginal problem), which is 
then determined algorithmically. This reduction applies more generally to symplectic geometry, 
relating invariant measures for a compact Lie group action to that for the maximal torus action; 
we state and prove our results in this more general symplectic setting. Our approach is in 
striking contrast to the existing solution to the computation of the supporting polytope by 
Klyachko and by Berenstein and Sjamaar, which made crucial use of non-Abelian features. 

In algebraic geometry, Duistermaat-Heckman measures correspond to the asymptotic distri- 
bution of multiplicities of irreducible representations in the associated coordinate ring. In the 
case of the one-body quantum marginal problem, these multiplicities include bounded height 
Kronecker and plethysm coefficients. A quantized version of the Abelianization procedure 
provides an efficient algorithm for their computation. 



The pure state of a quantum system is described by a vector in a complex Hilbert space, 
or more precisely by a point in the corresponding projective space. Herein we consider the 
finite dimensional case, for instance a spin system. Since the Hilbert space for multiple particles 
is given by the tensor product of the Hilbert spaces of the individual particles, its dimension 
grows exponentially with the number of particles. This exponential behavior is therefore the key 
obstruction to classical modeling of quantum systems. The observation is as old as quantum 
theory itself, and physicists ever since have tried to find ways around it. In that spirit, our paper 
presents an effective method to extract those physical features that "only depend on the one-body 
eigenvalues" associated to randomly-chosen quantum states of any fixed number of particles. Our 
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methods are rooted in geometry, and so further a nascent dialog between algebraic and symplectic 
geometry on the one hand, and the theory of quantum computation and quantum information 
on the other, with new results for each subject. We therefore begin with some context, and make 
an effort in the body of the paper to build a correspondence in terminologies. 

Typically, to address the aforementioned exponential complexity, physicists make use of a 
very simple yet powerful observation: Important properties such as energy and entropy often 
do not depend on the whole wavefunction but rather on only a small part, namely the reduced 
density matrix, or quantum marginal, of a few particles. For instance, the binding energy of 
a molecule is given as a minimization over two-electron reduced density matrices arising from 
A^-electron wavef unctions. Mathematically, the reduced density matrix is given as the contraction 
of (or trace over) the indices of the projection operator onto the wavefunction over the remaining 



particles (see (2.4) for the precise definition). The problem of characterizing the set of possible 
reduced density matrices, known as the quantum marginal problem in quantum information 
theory and as the A'^-representability problem in quantum chemistry |92[ [26| , has therefore been 
considered one of the most fundamental problems in quantum theory [98]. The general problem 
is computationally intractable, even on a quantum computer; more precisely it is QMA-complete 
and NP-hard |73| 174). However, the characterization of the one-body reduced density matrices of 
a pure global quantum state |6H [28| [62] admits a very elegant mathematical interpretation: it 
amounts to a description of the possible eigenvalues of the reduced density matrices; and this, 
just as for Horn's problem, requires the computation of moment polytopes for coadjoint orbits of 
unitary groups, as first observed in |23 l I28 | |6T]. An interesting consequence is that defining linear 
inequalities for the polytope can have physical interpretation: for instance, the Pauli principle is 
simply one linear inequality bounding the polytope of a fermionic system (see |25| [TT| [93l IMl |63] 
and |52| [53] [T3] for further examples). 

Random states play a fundamental role in physics. In classical statistical mechanics, the 
canonical state (or, "canonical ensemble") is the marginal probability distribution (Gibbs measure, 
or Boltzmann distribution) of states of the system arising from a uniformly random configuration 
of system and bath, subject to an energy constraint and fixed particle numbers for both system 
and bath [54i]. In quantum statistical mechanics, the canonical state of the system is the reduced 
density matrix of the uniform state on a subspace (encoding the energy constraint) of system and 
bath. In fact it can be shown that the canonical state almost always approximates the reduced 
density matrix of a random pure state in this subspace, as for large systems a concentration 
of measure occurs |89| 175} 135] . Earlier, and also out of thermodynamic considerations, Lloyd 
and Pagels in a seminal work computed the distribution of eigenvalues of the reduced density 
matrix of the system when system and bath are in a random pure state |76], i.e., in the simplest 
case of two large "particles", one representing the system and one the bath. In this paper we 
are concerned with computing exact eigenvalue distributions for any number of particles, rather 
than with asymptotic results or with simple model systems of two particles. More precisely, we 
are concerned with the following question: 

What is the probability distribution of the eigenvalues of the one-body reduced 
density matrices of a pure many-particle quantum state drawn at random from 
the unitarily invariant distribution? 

In this article, we answer this question completely by describing an explicit algorithm to compute 
the joint eigenvalue distribution of the reduced density matrices for an arbitrary number of 
particles of any statistics (distinguishable, Bose, Fermi). As a special case we easily recover the 
Lloyd-Pagels result for two distinguishable particles. As a corollary, our work naturally leads to 
a solution of the one-body quantum marginal problem in terms of a finite union of polyhedral 
chambers, thereby providing a complementary perspective on the work of Klyachko and of 
Berenstein-Sjamaar, and more recently of Ressayre, who each instead provided procedures based 
on geometric invariant theory (in particular the Hilbert-Mumford criterion) to list characterizing 
linear inequalities O [611 |9T] for the moment polytope. 

We emphasize that the computed eigenvalue distributions can be directly used to infer 
distributions of Renyi and von Neumann entropies, both of which play a fundamental role in 
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statistical physics and quantum information theory and are functions of the eigenvalues only. In 
particular, one can recover the average entropy of a subsystem [77' '86], which featured in an 
analysis of the black hole entropy paradox ^871 09] . For applications to the study of quantum 
entanglement, see 



52.5 



From a mathematical perspective, the eigenvalue distributions that we compute are Duister- 
maat-Heckman measures, which are defined using the push-forward of the Liouville measure on 
a symplectic manifold along the moment map (see §1.1 for the precise definition which we use in 



this article) [50| HSj HH |371 [38l EH]. The support of such a measure is a moment polytope, which 
in our physical context is the solution to the one-body quantum marginal problem. We will work 
and establish our results in this more general symplectic setting, in similar spirit to the existing 
solution to the one-body quantum marginal problem |9l El]. We remark that, in particular, the 
one-body quantum marginal problem subsumes the well-known Horn's problem, which asks for 
the possible eigenvalues of the sum A + B of two Hermitian matrices A, B with fixed eigenvalues 
[511 [601 l68l IM| [69] I70| . The corresponding eigenvalue distribution for randomly- chosen matrices 
A and B has also been studied in the literature, and our methods allow us to recover the main 



results of [2^ (Corollary 3.8); see also |32j for a more concrete approach 



However, one striking difference between existing approaches to the one-body quantum marginal 
problem and the corollary to our approach is that the subtleties associated with the non-Abelian 
nature of their solutions can be completely bypassed. For instance, their reliance upon cohomology 
of Schubert cycles, and the interplay of different Weyl groups and sub-tori that feature because 
of repeated use of the Hilbert-Mumford criterion, may be seen as incidental. The problem, 
even in the full generality of the symplectic setting, is at heart an Abelian one whose essential 
combinatorics is encoded in the maximal torus action together with taking finitely many explicit 
derivatives. Remarkably, this is more than philosophy, as it has real import for computation. 



The first mathematical contribution of this paper is an effective technique for providing explicit 
Duistermaat-Heckman measures under rather weaker assumptions than appear in the literature 



(§4). The fact that this subsumes our main question above, and in particular the one-body 



quantum marginal problem, while the existing literature does not, is the crucial second point 



of this work (§2.3). A third feature is our statement of the Abelianization procedure via the 



derivative principle for invariant measures (Theorem 3.3); though we have not seen this principle 



so formulated in the literature, there certainly are predecessors and the result should be equivalent 
to one of Harish- Chandra |45| . That a "quantized" version of our algorithm can be used to 



compute (efficiently, unlike existing algorithms, see §6.3) multiplicities for the branching problem 
of representation theory may be seen as the fourth mathematical consequence of our approach. 

The basic strategy we follow to address our main question above is a sequence of reductions. 
Firstly, our general quantum problem is replaced by an equivalent but more tractable one 



by "purification" of the quantum state (§2). This reduced problem is seen to satisfy a weak 



non-degeneracy assumption, so that the image of the moment map does not lie entirely on a wall 
in the relevant Weyl chamber. Under this assumption we can reduce via a derivative principle to 



the Duistermaat-Heckman measure for the maximal torus action (§3), which we evaluate by a 
"single-summand" algorithm along the lines of Boysal-Vergne [12J ( § 



14.2) 



This derivative principle holds for more general i('-invariant measures on the dual of the Lie 
algebra, t*: Every invariant measure can be reconstructed from its projection onto the dual of 
the Lie algebra of the maximal torus by taking partial derivatives in the direction of negative 
roots (Theorem 3.3). Again, we remark that such a reduction is not possible on the level of the 



supports of the Duistermaat-Heckman measures, i.e., on the moment polytopes. We also remark 
that in the case of the one-body quantum marginal problem, the Duistermaat-Heckman measure 
for the maximal torus action is the joint distribution of the diagonal entries of the one-body 
reduced density matrices (as opposed to their eigenvalues). Amusingly, this distribution can be 



viewed as the solution to the quantitative version of a classical marginal problem (§2.4). 



The single-summand algorithm alluded to above is an effective method for computing Duister- 



maat-Heckman measures for torus actions on projective spaces (Algorithm 4.17). Its name 
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stems from the fact that it amounts to evaluating a single summand of the kind that occurs 
in the well-known Heckman formula of Guillemin-Lerman-Sternberg [37], which expresses the 
measure as an alternating sum of iterated convolutions of Heaviside measures. We also describe 



an algorithm based on this latter formula (Algorithm 4.4) that can in particular be applied to 



projections of coadjoint orbits (Remark 4.6), and hence to the setting of Berenstein-Sjamaar [9j. 

Whenever a Hamiltonian group action can be quantized in a certain technical sense, Duister- 
maat-Heckman measures have an interpretation as the asymptotic limit of associated representa- 



tion-theoretic quantities |50| SSj |95l EHj |80l 1102) . In the second part of this article (§6 and §7), 
we thus study the representation theory connected to the one-body quantum marginal problem; 
here, the relevant multiplicities include the Kronecker coefficients, which play a major role in 
the representation theory of the unitary and symmetric groups [33], as well as in Mulmuley and 
Sohoni's geometric complexity theory approach to the P vs. NP problem in computer science 
|83[ IBH ESI [T^. It has been observed that the existence of a pure tripartite quantum state with 
given marginal eigenvalue spectra is equivalent to the asymptotic non- vanishing of an associated 
sequence of Kronecker coefficients |23| IGT] [22] , see also |28[ IFfj [TBI ITS] . For a similar connection 
in the context of Horn's problem and Littlewood-Richardson coefficients, see |72| EU] . 

Our main results in this context are quantized versions of our earlier theorems: Kronecker 
coefficients can be computed by applying finite difference operators to weight multiplicities which 
are related to the classical marginal problem: Instead of measuring the volume of a polytope, 
one has to count the number of lattice points in the polytope. This can be computed efficiently 
using Barvinok's algorithm [5j, and so leads to an efficient algorithm for computing Kronecker 
coefficients for Young diagrams with bounded height. Again, we shall establish the results in 
greater generality and recover the version for Kronecker coefficients as a special case. 

1.1. Notation and Conventions. Throughout this article, K will denote a compact, connected 
Lie group with maximal torus T <Z K, rank r = dimT, Weyl group W, respective Lie algebras t 
and t, and integral lattice A = ker exp |19tl57). We write ttk,t : ^* — )• t* for the projection dual to 
the inclusion t C t. We will think of weights as elements of the dual lattice. A* = Homz(A, Z) C t*, 
and identify a character x : T — t- U(1) with the weight dx/'^'^i G t*. We denote by d\ the Lebesgue 
measure on t* that is normalized in such a way that any fundamental domain of the weight 
lattice A* has unit measure. Let us also choose a positive Weyl chamber i'^ C t*; this determines 
a set of positive roots {a > 0} = {ai, . . . , or} Q t*. The set of negative roots is by definition 
{—a : a > 0}. We write tZ^Q for the interior of the positive Weyl chamber, which contains the 
strictly dominant weights. We will often identify t and its dual 6*, as well as t and t*, via some 
fixed X-invariant inner product ( — , — ) on t. 

For the special unitary group SU(d), which we always take to be the group of unitary d x d- 
matrices with unit determinant, we use the maximal torus consisting of diagonal matrices, on 
which the Weyl group Sd acts by permuting diagonal entries. The Lie algebra 5u{d) consists of 
anti-Hermitian matrices with trace zero, and our choice of invariant inner product is {X, Y) = 
— tr (XY) . Using it to identify t and t* , a positive Weyl chamber is given by the set of 
diagonal matrices A = diag(Ai, . . . , A^) with purely imaginary entries, summing to zero and 
arranged in such a way that iXi > . . . > iXd- This corresponds to choosing the positive roots 
o^j,ki^) = ii^j — A/fc) with j < k. The points in the interior i^Q are those A G with all distinct 
eigenvalues, i.e., iXi > . . . > iXd- 



Let M be a compact, connected Hamiltonian iiT-manifold of dimension 2n, with symplectic form 
UM and a choice of moment map : M — )■ t* [HI E]. The intersection Ak{M) = ^k{M) n t+ 
of its image with the positive Weyl chamber is a compact convex polytope, called the moment 
polytope or Kirwan polytope |42| l59l 140) . If M is a coadjoint iiT-orbit 0\ through some A G t!^, 
it will always be equipped with the Kirillov-Kostant-Souriau symplectic form and the moment 
map induced by the inclusion 0\ C 6*. Evidently, Ak{Ox) = {A}. 

Throughout this article, we will always impose the following non-degeneracy condition: 

Assumption 1.1. The moment polytope Ak{M) has non-empty intersection with the interior 
of the positive Weyl chamber, i*Q. 
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In view of |7H Lemma 3.9] and well-known facts about compact Lie group actions, this 
assumption in fact implies the following: The set ^~^{K ■ tt,Q) is an open, dense subset of M 
whose complement has Liouville measure zero. 



We show in §2.3 that Assumption 1.1 does not restrict the applicability of our techniques to 



the problem of computing eigenvalue distributions of reduced density matrices. 

The Duistermaat-Heckman measure DH^ is then defined as follows |30) : Push forward the 
Liouville measure = '^m/((2^)"'^0 ^ along the moment map compose with the 
push-forward along the quotient map Tx-i*^ i\ which sends all points in a coadjoint orbit 
Ox to A in the positive Weyl chamber; then divide the resulting measure by the polynomial 
PkW = na>o where p is half the sum of positive roots. That is, 

(1.1) DHiS = —{rM^Mm)- 

PK 



Assumption 1.1 ensures that DHf^ is a locally finite measure on the interior of the positive Weyl 
chamber. Its support is equal to the moment polytope. Note that pxW is equal to the Liouville 
volume of a maximal- dimensional coadjoint orbit 0\ |10| Proposition 7.26]. Therefore, the above 
is a natural definition to use in our context: It is normalized so that the Duistermaat-Heck- 
man measure associated with the action of on a generic coadjoint orbit 0\ is a probability 
distribution concentrated at the point A. 

If H is another compact, connected Lie group, with Lie algebra f), acting on M via a group 
homomorphism ip: H ^ K, then this action is also Hamiltonian, with a moment map given by 
the composition 

(1.2) = {dip)* o^k: M ^[)*. 

This in turn determines a Duistermaat-Heckman measure DH^. In particular, we can associate 
moment maps and Duistermaat-Heckman measures with all closed subgroups of K. In the case 
of the maximal torus T C if, we shall call the Abelian moment map and DH^ the Ahelian 
Duistermaat-Heckman measure^ in order to distinguish them from the non-Abelian moment map 
$x and the non-Abelian Duistermaat-Heckman measure DH^, respectively. Explicitly, 

(1.3) brIj = {<^T)4m) 

Throughout this paper we shall assume for simplicity that the Abelian moment polytope 
At{M) is of maximal dimension. This can always be arranged for by replacing T by the quotient 
meM where Tm denotes the T-stabilizer o f a point m g M . If T is the maximal torus 



of a semisimple Lie group, it follows already from Assumption 1.1 that Aj'(M) is of maximal 
dimension. As a consequence, T acts locally freely on a dense, open subset whose complement 
has Liouville measure zero \30\ Lemma 3.1]. In particular, generic points in M are regular for 
the Abelian moment map. Therefore the Abelian Duistermaat-Heckman measure is absolutely 
continuous with respect to Lebesgue measure on t*. By the Duistermaat-Heckman Theorem, 
DH^j in fact has a polynomial density function of degree at most n — r on each connected 
component of the set of regular values \30\ Corollary 3.3]. We shall call these components the 
regular chambers, and one can show that, except for the unbounded one, every such chamber 
is an open convex polytope. If the closures of two regular chambers have a common boundary 
of maximal dimension (i.e., of codimension one) then we shall say that the two chambers are 
adjacent and call the common boundary a singular wall. 

All Hilbert spaces which we consider in this article are complex and finite-dimensional. We 
write for the orthogonal projection onto a one-dimensional subspace Cip, and \\X\\2 = 
y/tT {X*X) = Y,. for the Hilbert-Schmidt norm of an operator X with singular values {sj). 
We use (— , — ) to denote inner products as well as the pairing between measures (or more general 
distributions) and test functions. We write 5p for the Dirac measure at p, i.e., the probability 
measure concentrated at the point p, and H^^ for the Heaviside measure which is defined by 
{H^, f) = /q°° f{tuj)dt. We sometimes use the letter P for probability distributions. 
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Throughout the paper when we speak of the quantum marginal problem we always refer to its 
one-body version as described in[§2] 

Finally, we offer a word of caution for people acquainted with the theory of geometric 
quantization |41 | I103j : Our quantum states do not arise via some quantization procedure from 
a classical symplectic phase space. In contrast, herein, as detailed in §2 below, the spaces of 
quantum states themselves are Hamiltonian manifolds. Probability distributions can be realized 
as quantum states of a special form, and the passage from quantum to classical is related 
to passing from a non-Abelian group to its maximal torus (see §2.4 for precise statements). 
The "semiclassical limit" well-known in geometric quantization does not have an analogous 
physical meaning in our setting; its significance is solely to connect the symplectic geometry with 
representation theory (see 



!6.1). 



2. Density Matrices and Purification 

The applicability of symplectic geometry to the quantum marginal problem relies on the close 
relation between the Lie algebra of SU((i) and the density matrices of quantum mechanics, and 
on the fact that restricting to certain subgroups has the physical meaning of passing to reduced 
density matrices, which describe the quantum state of subsystems. In this section we will describe 



this relationship in some detail (§2.1 §2.2), and show how one can reduce the general problem of 



computing joint eigenvalue distributions of reduced density matrices to the case of globally pure 



quantum states, that is, to the Duistermaat-Heckman measure for a projective space (§2.3). We 



briefly discuss how probability distributions and the classical marginal problem are embedded in 
our setup 



2.4) and describe some immediate physical applications (§2.5). 



2.1. Density Matrices. 

Definition 2.1. A density matrix is a positive-semidefinite Hermitian operator p of trace one 
acting on a finite-dimensional Hilbert space Ti. We will often choose coordinates and think of p 
as a matrix. If p is the orthogonal projection onto a one-dimensional subspace then we say that 
/9 is a pure state; otherwise, it is a mixed state. An observable is an arbitrary Hermitian operator 
acting on Ti. 

Density matrices on Ti describe the state of a quantum system modeled by the Hilbert space 
H: According to the postulates of quantum mechanics, the expectation value of an observable O 
is given by the pairing tr (Op) G M. Of course, p is characterized by these expectation values, 
even if we use anti-Hermitian observables instead and restrict to trace zero (since the trace of p 
is fixed). That is, we have an injection 

(2.1) p^{X^itv{Xp))£5u{-Hy 

which extends to an isomorphism between the affine space of trace-one Hermitian operators on 
Ti and 5u{T-L)* . This isomorphism is SU(?^)-equivariant; its inverse sends a coadjoint orbit 0\ to 
the set of Hermitian operators with eigenvalues 

(2.2) x, = ^r^+iX, (j = l,...,dim^), 

where the \j are the eigenvalues of A (eigenvalues are repeated according to their multiplicity). 
Let us choose coordinates Ji = and identify SU('H) = SU(d) accordingly. Then the eigenvalues 
Xi are just the diagonal entries of the matrix A G labeling the coadjoint orbit (see 
conventions) 



il.l 



for our 

It follows that (|2.2|) defines a bijection between the positive Weyl chamber tl and 



the set of eigenvalue spectra of trace-one Hermitian operators, which we think of as elements of 
the set {A G R'^ : Ai > . . . > A^, Y.j = 1}- 

Note that the set of pure states is identified with the coadjoint orbit through the highest 
weight of the defining representation of SU('H), that is, with projective space ¥{'H): The density 
matrix corresponding to a point [tp] G P('H) is simply the orthogonal projection onto Cip. 

Moreover, the Liouville measure on a coadjoint orbit Ox is identified via (2.1) with the unique 



SU(^)-invariant measure on the set of Hermitian matrices with spectrum A, normalized to total 
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Figure 1 . A Bloch sphere, its height function and the induced measure: a generic 
coadjoint SU(2)-orbit, its Abelian moment map and the Abehan Duistermaat- 
Heckman measure. 
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The non- Abehan Duistermaat-Heckman measure of Ox corresponds to the Dirac measure at A, 
while the Abelian Duistermaat-Heckman measure corresponds to the distribution of the diagonal 
entries of a density matrix with spectrum A chosen according to the invariant measure. 



Example 2.2. The Lie algebra 5u(2) is three-dimensional, generated by 
ax, CTy, Gz- Functionals in its dual su(2)* can be identified with points in 



2 times the Pauli matrices 

3 



by evaluating them 

at these generators. In this picture, the inverse of (2.1) associates to a vector r = (x,y,z) the 
Hermitian matrix p{r) = ^(1 + r • <?), where (J is tlie PcLuli vector ay,az)- The z- axis is 
identified with the Lie algebra of the maximal torus, Mctz, its positive half-axis with our choice of 
positive Weyl chamber t^^, and (0, 0, 2) with the corresponding positive root a > 0. The coadjoint 
action of elements in SU(2) amounts to rotating the Bloch vector r via the two-fold covering 
map SU(2) — )• S0(3). Therefore, coadjoint orbits are spheres, commonly called Bloch spheres in 
quantum mechanics. They can be labeled by their radius r, that is, by their intersection with 
the positive half of the z-axis. Points on such a sphere correspond to Hermitian matrices with 
eigenvalue spectrum i-^-^, ^-^)- Note that p{r) is positive-semidefinite (i.e., a density matrix) if 
and only if r is contained in the unit ball of M.^. 

The non- Abelian moment map is just the inclusion map of a Bloch sphere into M"^. Hence its 
composition with the quotient map Tsu(2) sends all points in a Bloch sphere of radius r > to 
r, while the Abelian moment map projects all points onto the z-axis. The Liouville measure is 
equal to the usual round measure, normalized to total volume r. Therefore, the non-Abelian 
Duistermaat-Heckman measure of a Bloch sphere with radius r is equal to the Dirac measure 6r, 
while the Abelian Duistermaat-Heckman measure is obtained by pushing forward the Liouville 



measure onto the z-axis (see Figure 1 ). As already observed by Archimedes, any two zones of 
the same height on a sphere have the same area. Hence this latter measure is proportional to 
Lebesgue measure on the interval [— r, r]. An analogous statement holds for arbitrary projective 



spaces (Proposition 4.10) 



Remark 2.3. Observe that the components (X, $gu(-H)) ■ C'a — s- M of the moment map send a 
quantum state to the expectation value of the corresponding observable —iX. Without loss 
of generality, we may assume that X generates a one- dimensional torus and that X has one- 
dimensional eigenspaces. Then {X, $gu(-^)) is just the moment map for the action of the torus 
generated by X and its distribution can be computed immediately by using the Abelian Heckman 
formula (Theorem 4.1). This gives a short and conceptual proof of the formula derived recently 
in [99J. 



2.2. Reduced Density Matrices. Composite quantum systems are modeled by the tensor 
product of the Hilbert spaces describing their constituents. It is useful to think of these subsystems 
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as individual particles, although they can be of more general nature; for instance, the subsystems 
can describe different degrees of freedom such as position and spin. Depending on whether the 
particles are in principle distinguishable or indistinguishable, we distinguish two basic classes of 
composite systems, which are of fundamentally different nature. 

If the quantum system is composed of N distinguishable particles, its global quantum state 
is described by a density matrix on the tensor product Ti = Tii ■ ■ ■ <E> T-Ln, where the Tij 
are the Hilbert spaces describing the individual particles. Quantum mechanics also tells us 
that observables Oj acting on a single subsystem Tij correspond to tensor product observables 
l'Hi(Si...®'Hj_i ^ Oj (8) l-Hj+i®...®'Hjv) which act by the identity on all other subsystems. By non- 
degeneracy of the inner product, there exists a unique density matrix pj on Tij such that 

(2-4) tr Oj (g) 1hj+^®...®Hn)p) = (OjPj) 

for all observables Oj. It describes the quantum state of the j-th subsystem. 

Definition 2.4. The density matrix pj is called the (one-body) reduced density matrix or quantum 
marginal for the j-th particle of the quantum system. 

Note that we can embed S\]{'Hj) into SU(^) by Uj i-)- 'i-ni(S)...(S)'Hj-i ® Uj (8) l-Hj+i^.-.^Wiv This 
induces an embedding on the level of Lie algebras. The dual projection su('H)* — )• su('Hj)*, given 



by restricting functionals to the subalgebra, is identified by (2.1) with the map /9 1— )• pj. Similarly, 
the group homomorphism from the Cartesian product SU(?^i) (8) ... (8) S\]{7iN) to S\J{T-L) given 
by (C/i, . . . , Un) I— 7- C/i (8 . • • (8 Un induces the map p i— )• (pi, . . . , /Otv) sending a density matrix to 
the tuple of all its one-body reduced density matrices. 

The (one-body) quantum marginal problem for distinguishable particles asks for the possible 
tuples of one-body reduced density matrices (pi, . . . , p^) of an arbitrary density matrix p with 
fixed spectrum, or, equivalently, for the possible tuples of their eigenvalues. By the above 
discussion, this is precisely equivalent to determining the moment polytope Ak{M) associated 
with the Hamiltonian action of the subgroup K = SU(?^i) (8) ... (8" SU(^7v) on a coadjoint orbit 
M = for Sl](T-L), with moment map as defined in (1.2). The quantum marginal problem 



for globally pure states is the special case where M = ¥('H). Moreover, the joint eigenvalue 
distribution of reduced density matrices we set out to compute in this article corresponds to the 
non-Abelian Duistermaat-Heckman measure DH^ as defined in (1.1): Up to the identification 



A I— 7- A between spectra of trace-one Hermitian operators and the positive Weyl chamber as 



defined in (2.2), it is given by 



(2.5) {tkW^k)* ( = -4yjPk DHi^ 



volMJ volM 

Note that we divide by the Liouville volume of M, which is just Ps\JCH)W^ ^o obtain a probability 
measure. 

Similarly, the joint distribution of the diagonal entries of the reduced density matrices corre- 



sponds to the Abelian Duistermaat-Heckman measure DHJ 



If the quantum system is composed of indistinguishable particles, each particle is of course 
modeled by the same Hilbert space Tii. The global state of the system is described by a density 
matrix supported on an irreducible sub-representation 7i C T-lf^, namely 7i = Sym^('Hi) 
for bosons and ((Hi) for fermions (but we can in principle also consider other irreducible 
sub-representations which correspond to more exotic statistics). Note that since every such 
density matrix commutes with permutations, all the one-body reduced density matrices are equal. 

Note that we can let single-particle observables O act more intrinsically by the symmetric 
expression jj{0 (8> l^®(]v-i) -|- . . . -|- l^®(jv-i) ® O), without changing their expectation values. 

Up to a factor N , this corresponds to the embedding of Lie algebras induced by the diagonal 
map U ^ U ® . . . ®U , which is of course precisely the action of SU(?^i) on the representation 
H. This embedding therefore induces the map p i— )■ pj = Npi in the same way as described 
above. It follows that the (one-body) quantum marginal problem for indistinguishable particles 
amounts to determining -^A/f (M) for the induced action of K = SU(^i) on a coadjoint orbit 
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Setting 


Hilbert space Ti 


Group K 


N distinguishable particles 


. . . (g) C"*^ 


SU(di) X ... X SU(djv) 


bosons 




SV{d) 


N fermions 


A^(C'^) 


SU((i) 



Table 1 . The quantum marginal problem is modeled by the action of the group 
on a coadjoint SU(?^)-orbit M = Oy 



M = of S\J{T-L), and that, up to the identification A i— )• A, the eigenvalue distribution of the 
reduced density matrix is given by 

(2.6) K,{TKU'^>Kh f-^) = -^PK^*mM), 



volMJ volM 

where the linear map k{X) = counteracts the factor in the moment map. 

Remark 2.5. In |Table~T| we have summarized the groups and spaces relevant for the quantum 
marginal problems of main physical interest, and we will focus on these in the remainder of this 
article. One can also combine both cases, e.g., to describe a quantum system composed of two 
different sorts of indistinguishable particles (as happens for the purified double of a bosonic or 



fermionic quantum marginal problem as defined in §2.3 below), or a number of indistinguishable 
particles each of which have multiple internal degrees of freedom. In the latter case, arbitrary 
irreducible representations of the special unitary group can appear if one restricts to the reduced 
density matrices corresponding to only some of the degrees of freedom (see, e.g., |64] ) . 

2.3. Purification. Let T-L be an arbitrary finite-dimensional Hilbert space. It is well-known that 
every density matrix /? on ^ is the reduced density matrix of a pure state in [ip] G ¥(T-L (8) Ti), 
called a purification of the quantum state p. Indeed, if p = ^^PiPvi is the spectral decomposition 
of p then we can simply choose ip = y/PiVi Vi. In this sense, the global state of a quantum 
system can always be described by a pure state; reduced density matrices occur only in the 
description of the states of its subsystems. This motivates the following definition: 

Definition 2.6. For any unitary ^-representation Ti, we define the purified double to be the 
Hamiltonian K x i^-manifold ¥('H (S" T-L), where K = SU(^), equipped with the moment map 
constructed in the usual way by embedding into u('H C3 T-L)* and "restricting" the functionals to 
elements in t(Bt- 



Observe that if Ti is one of the representations of |§2.2 modeling a setup of the quantum 



marginal problem, then the purified double corresponds to the pure-state quantum marginal 
problem where one has adjoined a single distinguishable particle modeled by 7i. 

The purification [ip] G ¥(T-L (g) T-L) of a quantum state p on 7^ is unique up to a unitary acting 
on the second copy of H. Evidently, such operations do not change the reduced density matrix 
p = (-P,/,)i and they leave the eigenvalue spectrum of (-P^)2 invariant. In particular, the eigenvalue 
spectra of the reduced density matrices (-P^)i and {P^)2 are always equal. This implies that 
we can reduce the quantum marginal problem to the case of globally pure states, both for 
distinguishable and indistinguishable particles: 

Proposition 2.7. Let be a coadjoint orbit of K = SIJ('H), with X i*^ corresponding 
to the eigenvalue spectrum of a density operator. Then A G A/^(0^) if and only if (A, A) G 

In other words, 

A^^j^im^-H)) = U Ak{0-,) X {A}, 
AeA 

where A = {A G t!^ : Xj > 0} is the convex subset of the positive Weyl chamber corresponding to 
the eigenvalue spectra of density operators. We can similarly reduce the problem of determining 



10 



M. CHRISTANDL, B. DORAN, S. KOUSIDIS, AND M. WALTER 



the joint eigenvalue distribution to the case of globahy pure states: For this, let us define 
probability measures 



(2.7) 



MP(-H(g)-H) 
KxK)*\'^KxKl* \ volP(^®7^) 



voIO^ 



where vol denotes the Liouville volume. 



Proposition 2.8. The measures in (|2.7|) are related by 

(P,/) 



1 
Z 



dXpj^iX) {P-„f{-,X)), 



for all test functions f € Cf,(t!j_ ® t^J,), where dX is Lebesgue measure on and Z a suitable 
normalization constant. 

Proof Each of the one-body reduced density matrices of a Liouville-distributed bipartite pure 
state in ¥(T-L 71) is distributed according to the Hilbert-Schmidt measure restricted to the set 
of density matrices. In particular, its eigenvalues are distributed according to the well-known 
formula of [76l 1105] , so that 



^p\(X) l^{X)d\ 



with 1^ the indicator function of A and Z a suitable normalization constant. We have just seen 
that both reduced density matrices necessarily have equal eigenvalue spectrum. This implies that 



KxK!* 



volP(^«)^) 



,9) 



1 



dX 



9- 



See Corollary 5.17 for an independent derivation using the techniques of this paper. It follows 
that 



dX 



1 
z 
1 
z 
1 
z 



J) 



dX pf,{X) 



{TK<fKr{f{-,X) 



dXp\{X) {P-^,fi-,~X)). 



□ 



Note that P^^, and in particular the eigenvalue distributions (2.5) and (2.6), vary continuously 
with the global spectrum A. Proposition 2.8 therefore implies that we can reconstruct them from 
the eigenvalue distribution for the purified double by taking limits. If the latter distribution has 
a continuous Lebesgue density, as will often be the case, then we can simply restrict this density 
function to the global spectrum A of interest. 



We will now show that Assumption 1.1 is always satisfied when working with the purified 
double. In quantum-mechanical terms, we have to show that there exists a global pure state 
[ip] E ¥(T-L 0T-L) such that the eigenvalue spectra of all the reduced density matrices are non- 
degenerate (with respect to the quantum marginal problem where we have added a single 
distinguishable particle with Hilbert space H). 



EIGENVALUE DISTRIBUTIONS OF REDUCED DENSITY MATRICES 



11 



For distinguishable particles, where T-L = C^^ (g) . . . (8) C"^^ , this follows from the following 
general criterion, since the purified double is constructed by adding an additional Hilbert space 
of dimension d7v+i = dimTi = di - ■ ■ d^: 

Lemma 2.9. Let N > 1 and di < ... < djsf < cZtv+i- Then there exists a global pure state 
in P(C'^^ (8) ... (8) C'^'^ (8 C^^+i) whose one-body reduced density matrices have non-degenerate 
eigenvalue spectra if and only if 



+ 1. 



Proof. The condition is clearly necessary, since it follows from the singular value decomposition 
that at most YiiLi di eigenvalues of pN+i can be non-zero. 

For sufficiency, let us construct a state with the desired property: For this, we consider the 
standard tensor product basis vectors (8 . . . (8 Cj^ of C^^ (8 • • • <8 C^^, labelled by integers 
ij £ {1, . . . ,dj}, j = 1,...,N. We choose a subset of (iAr+i — 1 many such basis vectors 
Cii(fc) 8> . . . 8) &i!^{k) ill such a way that, for each subsystem j = 1, . . . , N , at least dj — 1 of the dj 
integers occur. This clearly is possible by our assumptions. Finally, we set 

cfjv+i — 1 

V'= XI 2"''eij(A.) (8 . . . (8 ei^(fc) (8 Cfc. 
ifc=l 

Then [ip] is a pure state such that all of its one-body reduced density matrices have non-degenerate 
eigenvalue spectrum. □ 



For bosons and fermions, we will use the following lemma to show that [Assumption 1.1 is 
satisfied: 

Lemma 2.10. The convex hull of the weights o/Sym^(C'^) has maximal dimension. The same 
is true for (C^) if N < d. 

Proof. In the following we write ujj for the weight corresponding to the character diag(ti, . . . , t^) 
tj , with j = 1, . . . ,d. 

(1) Sym^(C'^): The vectors Cj 8> • ■ • 8> ej are weight vectors of weight Nooj, with j = 1, . . . ,d. 
Clearly, the convex hull of these weights already has maximal dimension. 

(2) A^(C'^), where 1 < < d: It is well-known that the weights are given by "^j^jOJj for 
all A^-element subsets J C {1, . . . ,d}. (The corresponding weight vectors are the well-known 
occupation number basis vectors for fermions.) Fix any such weight, say, the one corresponding 
to J = {1, . . . , A^}. The difference vectors between this weight and the weights obtained by 
replacing a single element of J are proportional to the positive roots aj^kW = ii^j — ^k) with 
j £ J and k G {1, . . . ,d}\J. There are at least d — 1 such roots, and they form a basis of t*. □ 



Lemma 2.10 implies that [Assumption 1.1 is also satisfied for the purified double of the bosonic 



and fermionic quantum marginal problems: Indeed, it clearly suffices to show that in each case 
there exists a density operator p onH such that both p and its one-body reduced density matrix 
pi have non-degenerate eigenvalue spectrum (then any purification of p has the desired properties). 



By [Lemma 2.10 there exists a convex combination of weights YlkPk^k £ t>o- By perturbing 
slightly, we can arrange for the weights (pk) to be mutually disjoint. Choose corresponding 
(orthogonal) weight vectors Vk £ Ti and consider the density matrix p = YlkP^^Vk- Clearly, both 
p and its one-body reduced density matrix pihave non-degenerate eigenvalue spectrum. 

To summarize, we have shown that the problem of computing the joint eigenvalue distribution 
of reduced density matrices is equivalent to the computation of Duistermaat-Heckman measures 



associated with certain Hamiltonian group actions (cf. Table 1). Moreover, by passing to the 
purified double, we can always reduce to the case where M = F(T-L) is a projective space satisfying 
Assumption 1.1 
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2.4. Probability Distributions. Under the identification (2.1), elements of the dual of the Lie 
algebra of the maximal torus correspond to diagonal density matrices. These are precisely the 
diagonal matrices with non- negative entries summing to one, and can therefore be interpreted as 
probability distributions of a random variable Z with values in the orthonormal basis (ej) we have 
chosen. This interpretation is in agreement with quantum mechanics: If we perform an actual 
measurement of a density matrix p with respect to this orthonormal basis then the probability of 
getting outcome is given precisely by the diagonal element tr (Pe-p) = {ei,pei) = pi^i. 

Note that the moment map for the action of the maximal torus T C S\J{'H) on the projective 
space ¥{7i) corresponds to sending a pure state [■0] onto its diagonal. As we vary [■(/;] over all pure 
states in IP(^), the diagonal entries attain all possible probability distributions. In other words, 
the Abelian moment polytope Aj,{¥(T-L)) is just the simplex A defined in §2.3 The corresponding 



Duistermaat-Heckman measure is equal to a suitably normalized Lebesgue measure on A (this 



is a special case of Proposition 4.10 below). 



Now consider as in §2.2 the case of distinguishable particles. Choose orthonormal bases 
to identify Tik = C^*, and therefore ^ = (g) . . . ® C^^ using the tensor product basis. Note 
that we can interprete diagonal density matrices /9 on as the joint probability distribution of a 
tuple of random variables (Zi, . . . , Zj\[), where each takes values in the standard basis of the 
corresponding C'^'', by setting 

P(Zi = Cii, . ..,Zn = Ci^) = tr (Pei^C>5...®e,^P 

The marginal distributions of the random variables Z^ in the sense of probability theory are then 
given by 



^l,...,^k,■■■,^N 



where for the second identity we have used that p is a diagonal matrix. That is, the marginal 
distributions of the Z^ are precisely described by the reduced density matrices pk (i.e., by the 
quantum marginals), which are also diagonal if p is diagonal. 

Accordingly, the moment polytope Ar(P(?^)) for the action of the maximal torus T C 
SU((ii) X . . . SU((iiv) on the set of pure states describes the tuples of marginal probability 
distributions that arise from joint distributions of the (Zi, . . . , Z]\f). This (univariate) classical 
marginal problem is of course trivial, since there are no constraints on the joint distribution. 
However, its quantitative version, which corresponds to computing the Abelian Duistermaat- 
Heckman measure DHp^^^, is interesting and not at all trivial to solve. In fact, the problem 
of computing joint eigenvalue distributions of reduced density matrices, which we set out to 
solve in this article, can be reduced to the computation of DHp^^^^. This reduction, or rather the 
generalization which we describe in ^ below, is at the core of the algorithms presented in |4 



2.5. Physical Applications. As indicated in the introduction, the eigenvalue distributions 



(2.5) and (2.6) have direct applications to quantum physics. In quantum statistical mechanics, 
among others, one typically studies bipartite setups Ti = Tis ®%e composed of a system S 
and an environment (or bath) E. Randomly- chosen pure states give rise to a distribution of 
reduced density matrices p^, whose properties vary with the size of the environment. Physical 
motivations have lead to the computation of the corresponding eigenvalue distribution [76], which 



we can easily re-derive using the techniques of this paper (Corollary 5.15). Note that many basic 



physical quantities are functions of the eigenvalues, such as the von Neumann entropy 
H{S) = H{ps) = -tr {ps log ps) = -^^J log 
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where (Xj) are the eigenvalues of ps, or more general Renyi entropies and purities (cf. Corol- 



lary 5.13). The average von Neumann entropy of a subsystem |77| |86] in particular has featured 
in the analysis of the black hole entropy paradox [49j. 

We can also consider other coadjoint orbits such as Grassmannians: Here, the density matrix 
corresponding to a (i- dimensional subspace Ti' C Jij^ (^^He is the normalized projection operator 
p = 1-^/ /d, and the reduced density matrix p^ is interpreted as a canonical state in the sense 
of statistical mechanics |89l I75ll35] . The probability distributions we compute can therefore be 
used to analyze the typical behavior of canonical states. 

The tripartite case, in itself already interesting from the perspective of the quantum marginal 
problem, is also highly relevant to applications: It corresponds to the situation where S itself 
is composed of two particles A and B, so that H = Ha ®T~Lb®T~Le- In the study of quantum 
entanglement, remarkable recent progress has been made by analyzing the entanglement properties 
of the two-body reduced density matrix pAB of a randomly-chosen pure state in large dimensions, 
where the concentration of measure phenomenon occurs |47| HEl HI El EZI- In particular, a 
negative resolution of the additivity conjecture of quantum information theory |94| has recently 
been obtained by related methods [l6| [T]. The joint eigenvalue distribution of the reduced 
density matrices in particular determines quantum conditional entropies and quantum mutual 
informations, that is, the quantities 

H{A\B) = H{AB) - H{B) = H{E) - H{B), 

I{A : B) = H{A) + H{B) - H{AB) = H{A) + H{B) - H{E), 



since the eigenvalue spectra of pAB and pE are equal (cf. §2.3). They have immediate applications 
to entanglement theory; for example, the quantum mutual information provides an upper bound 
on the amount of entanglement that can be distilled from a quantum state |24| . 

In all these applications, most known results are for large Hilbert spaces, since the techniques 
employed rely on asymptotic features such as measure concentration. Our algorithms require 
no such assumption. In particular, they are well-suited for low-dimensional systems, which 
previously remained inaccessible. 



3. Derivative Principle for Invariant Measures 



In this section we will describe a fundamental property of ET-invariant measures on t* that 
are concentrated on the union of the maximal-dimensional coadjoint orbits (that is, on K ■ tl^g). 
Every such invariant measure can be reconstructed from its projection onto t* by taking partial 



derivatives in the direction of negative roots (Theorem 3.3). In particular, this implies that the 



non-Abelian Duistermaat-Heckman measure DH^ can be recovered from the Abelian Duister- 



maat-Heckman measure DH^^ ( [Corollary 3.4| . 

For the invariant probability measure supported on a single coadjoint orbit of maximal 
dimension, this follows from a well-known formula of Harish-Chandra, as was already observed 
by Heckman: 

Proposition 3.1 

Then, 



45] Theorem 2], [50l (6.5)]). Let Ox be a coadjoint orbit through A G tt^Q. 



Us- 



DH 



where the partial derivatives and the restriction are in the sense of distributions. 

Proof. Harish-Chandra's formula for the Fourier transform of a coadjoint orbit states that 



l{w) i(wX,X) 



n 

a>0 



1 

i{a, X) 
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for every X E t which is not orthogonal to a root (see [10, Corollary 7.25] for a recent account). 
Here, /(t/;) is the length of the Weyl group element w. This implies that the Abelian Duister- 
maat-Heckman measure is given by the following alternating sum of convolutions 



DH 



(3.1) 

wew 

Recall that H^^ is the Heaviside measure defined in 
fundamental theorem of calculus, we have du,H,„ = 6n. I'herefore 



il.l 



by (H^J) = f{tu:)dt. By the 



(3.2) 



DH 



\a>0 



and the assertion follows if we restrict to the interior of the positive Weyl chamber. □ 
Example 3.2. Every Bloch sphere of radius r > is a coadjoint orbit of maximal dimension (cf. 



Example 2.2). We have seen that DH^^ is equal to j,](z)(i2;, where dz is Lebesgue measure 

on the z-axis. In agreement with Proposition 3.1[ we observe that 



ti>0 



5r- 



dzl[-r,r]{z)dz\ 

Theorem 3.3. Let v he a K -invariant Radon measure on t* satisfying z^(t* \ K ■ t^g) = 0. Then, 

1 



Jl d-a {t^k,t)*{t^) 



\a>0 



PK 



'>0 



where the partial derivatives and the restriction are in the sense of distributions. 

Proof. Let / G C^(t>o) ^ test function, which we extend by zero to all of t*, and set 
9 '■= {'^k,t)* {{Ylayo^a) /)■ By definition and assumption, respectively. 




J) = {v^a) = {^\ 



9)- 



Since is a X-invariant measure, we can use Fubini's theorem to replace g by its ii'-average. On 
each maximal-dimensional coadjoint orbit 0\ C K ■ iZ^Q, this average is given by 



1 



volO 



1 



A 



(/i, 



Oa 




1 



\q>0 



which by Proposition 3.1 is precisely equal to /(A)/px(A). In other words, the averaged function 



is on X . tI,Q equal to the pullback (r^)* {f /pk)- We conclude that 

{t^lx-e '9) = W\K.i* ^(tk)* (—]) = {—irK)*{i^) 
'■^ So So \Pk J Pk 



,/>■ 



□ 



Corollary 3.4. The Duistermaat-Heckman measures as defined in §1.1 are related by 









\a>0 / 


t>o 





Proof. I Assumption 1.1 guarantees that we can apply |Theorem 3.3 to the push-forward of the 
Liouville measure along the non- Abelian moment map ^k- D 

This is the derivative principle alluded to in the title of this section. As we shall see in the 
following, it is a powerful tool for lifting results about the Duistermaat-Heckman measure for 
torus actions to general compact Lie group actions. 

Remark 3.5. According to |104t §3.5], [Corollary 3.4 was already known to Paradan and also 
follows from a different result of Harish-Chandra. In |6.2 we will describe another way to establish 
it by using the connection between Duistermaat-Heckman measures in algebraic geometry and 
multiplicities in group representations. 
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Remark 3.6. Note that Theorem 3.3 completely determines the measure z/ from its projection onto 
t*, since i' is by assumption concentrated on the union of the coadjoint orbits of maximal dimension. 
Similarly, the n on-Abelian Du istermaat-Heckman measure DH|^ can be fully reconstructed from 



DH^^ by using Corollary 3.4 



We stress that it is oftentimes not necessary to explicitely compute the non-Abelian Duister- 
maat-Heckman measure. Indeed, [Corollary 3.4 is of course by definition equivalent to 



n /) 



for all / E C^(t^g), so that we can reduce the computation of averages over DHf^ directly 



to integrations with respect to the Abelian Duistermaat-Heckman measure (cf. proof of Corol- 



lary 5.13) 



Remark 3.7. It follows from Corollary 3.4 and the discussion in §1.1 that, on each (open) regular 



chamber, the non-Abelian Duistermaat-Heckman measure also has a polynomial density, namely 
the partial derivative in the directions of the negative roots of the density of the Abelian measure. 
However, there could still be non-zero measure on the singular walls separating the regular 
chambers. If we would like to exclude this then we need to understand the smoothness properties 
of the Abelian density function in the vicinity of singular walls, or, equivalently, the nature 
of the term by which the polynomial density changes when crossing a singular wall. If this 
jump term vanishes to order at least R on the wall, then the Abelian density function is at 
least i?-times weakly differentiable in the vicinity of the wall, and therefore the non-Abelian 
Duistermaat-Heckman density is also absolutely continuous there. This vanishing condition can 
be checked explicitly for each singular wall using the jump formula described in |§4.2[ 

In case the vanishing condition is satisfied, the non-Abelian moment polytope Ak{M) is equal 
to the closure of a finite union of regular chambers for the Abelian moment map: Indeed, on 
each regular chamber the density polynomial is either equal to zero, or it is non-zero on an open, 
dense subset. 



We cannot resist giving an easy application of Corollary 3.4 to the generalized Horn's problem: 
Here, the goal is to describe the sum of two coadjoint orbits 0\ + (Horn's original problem 
referred to coadjoint orbits of the special unitary group). In other words, one considers the 
diagonal action of K on Ox x O^, which is Hamiltonian with moment map {X, Y) X + Y, and 
one would like to describe the associated moment polytope or Duistermaat-Heckman measure. 

Corollary 3.8 ([29j). Let A G and fi £ i*^. Then, 



DH 



K 



5](-1)'(-)5^a*DHS^, 



where l{w) is the length of the Weyl group element w. 
Proof. Clearly, since T is Abelian, 



DHg^^DHg^ 



Since X + fi £ tig. Assumption 1. 1 



assertion follows together with (3.2) 



is satisfied. Therefore, Corollary 3.4 is applicable, and the 



DE 



K 



■wX 



□ 



The general case where both A and fi are contained in the boundary of the positive Weyl 
chamber can be treated as in |29) by taking limits. Of course we can also expand DH^ as an 



alternating sum of convolutions by using (3.1) or its version for lower-dimensional coadjoint 
orbits [ini Theorem 7.24]. 
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4. Algorithms for Duistermaat-Heckman Measures 
In this section we present two algorithms for computing Duistermaat-Heckman measures. 



Both algorithms are based on the derivative principle from §3 , in that they first compute the 
Abelian measure and then take partial derivatives according to Corollary 3.4[ 

The first algorithm, the Heckman algorithm, is based on the Heckman formula by Guillemin, 
Lerman and Sternberg, which expresses the Abelian measure as an alternating sum of iterated 
convolutions of Heaviside measures. The density function of each such convolution is piecewise 
polynomial and can be evaluated inductively using recent work of Boysal and Vergne. While 
very useful for computing low-dimensional examples, the resulting algorithm is rather inefficient 
due to the large number of summands. 

Our second algorithm, the single- summand algorithm, is based on another formula for the 
Abelian Duistermaat-Heckman measure in the case where M is the projective space of an arbitrary 
finite-dimensional representation. It turns out that this formula is equivalent to evaluating a 
single iterated convolution of the above form (hence the name of the algorithm). It can therefore 
be computed in a similar way, but much more efficiently. Since by passing to the purified double 
the quantum marginal problem can always be reduced to the case where M is a projective space 



(§2.3), this solves the problem of computing eigenvalue distributions of reduced density matrices 



in complete generality. 

4.1. Heckman Algorithm. Before stating the Heckman formula by Guillemin, Lerman and 
Sternberg, let us recall the following renormalization process as described in |37]: 

Suppose that there are only finitely many fixed points of the action of the maximal torus T on 
M. For each such fixed point p G M^, consider the induced representation of T on the tangent 
space TpM. The weights of this representation are called isotropy weights and we can always 
choose a vector 7 G t* which is non-orthogonal to all isotropy weights (for all tangent spaces). 
The process of multiplying by —1 those isotropy weights that have negative inner product with 7 
is then called renormalization, and the resulting weights are called renormalized weights. See 



Remark 4.5 for a discussion of the case where M is a projective space and §5 for examples. 

Theorem 4.1 (|37|). Suppose that there are only finitely many torus fixed points p G M'^ . Denote 
by Up the number of isotropy weights in TpM that are multiplied by —1 during renormalization 
and by ujp^i, . . . ,u]p^n the resulting renormalized weights. Then, 

DRl = i-^r^S'S'TiP) * 

p€MT 

with Hi2i the Heaviside measure defined by {H^, f) = dtfiujt). 

In other words, the stationary phase approximation for the Fourier transform of an Abelian 
Duistermaat-Heckman measure is exact. This generalizes the Harish-Chandra formula for 



coadjoint orbits (3.1), which we used to establish Proposition 3.1 

Observe that each summand of the Heckman formula can be written as the push-forward of 
the standard Lebesgue measure dt on ]R>q along a linear map of the form P : [tj.) 1— )• X]fc=i ^fc'^fe) 
translated by ^t{p)i since 

(4.1) H^,*...*H^^=P4He,*...*HeJ = P4dtL„ ). 

'"■>o 

In a recent paper [12], Boysal and Vergne have analyzed general push-forward measures of this 
form under the assumption that the vectors ujk span a proper convex cone (i.e., a convex cone of 
maximal dimension that does not contain any straight line). This ensures that the measure is 
locally finite and absolutely continuous with respect to Lebesgue measure on t* . This assumption 
is certainly satisfied for the renormalized isotropy weights occurring in the Heckman formula (by 
the very definition of renormalization and our assumption that the Abelian moment polytope 
has maximal dimension). 

Let us briefly review their results: It is well-known that the push-forward measure has a 
piecewise homogeneous polynomial density function of degree n — r. Here, the chambers are 
the connected components of the complement of the cones spanned by at most r — 1 of the 
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weights (wfc). Except for the unbounded chamber, they are open convex cones. Walls are by 
definition the convex cones spanned by r — 1 hnearly independent weights Similarly to §1.1 if 



the common boundary of the closure of two chambers is of maximal dimension then this common 
boundary is a wall; moreover, every wall arises in this way. Note that the union of the walls is 
precisely the complement of the union of the chambers 

Let A± be two adjacent chambers which are separated by a wall W, and choose a normal 
vector ^ E t* pointing from A_ to A+. Order the weights such that precisely cji, . . . , Um lie on 
the linear hyperplane spanned by W. In the following, we shall freely identify differential forms 
and the measures induced by them. Denote by dw the Lebesgue measure on the hyperplanes 
parallel to W, normalized in such a way that 

(4.2) dX = dwA di 

where is the pullback of the standard volume form of M along the coordinate function 
{—,(,)■ Denote by f± the homogeneous polynomials describing the density function / on 
A-I-. Finally, consider the push- forward of Lebesgue measure on M>q along the linear map 
^w'- (uk) I— ^ X^fcLi ^yf^fc- Its density with respect to dw is given by a single homogeneous 
polynomial on the wall W, since W is always contained in the closure of a chamber for P^. 
Denote by f^^, any polynomial function extending it to all of t*. Then the result of Boysal and 
Vergne is the following | ,12, Theorem 1.1]: The jump of the density function across the wall is 
given by 

where Res |^_q5 = Q-i is the residue of a formal Laurent series g = "^^k^k^^- (The residue 
appears as part of an inversion formula for the Laplace transform.) 

In the case where only a minimal number of weights lie on the linear hyperplane spanned by 
W (m = r — 1), the wall polynomial /^^ can be chosen as a constant, since the corresponding 
push-forward map is merely a change of coordinates: 

Lemma 4.2. Suppose that precisely r — 1 weights wi, . . . ,u)r-i He on spanVl^. Then, 



Proof. Since the map P^j^ : — )• span W^, (uk) ^ X]fc=i ^fc'^fc along which we push forward 
is a linear isomorphism, the polynomial can be chosen as the constant of proportionality 
between the push-forward of Lebesgue measure on W'^'^ and the measure dw. We can compute 
its value by comparing the volume of the parallelotope spanned by the ((2;^) with respect to the 
two measure. For the former measure, this is of course one, while for the latter it follows from 



( |4.2D that 

dXiiUi, . . .,idr-i,i) = dw{ibi, . . . , Wr-i) IkiP- □ 

This immediately gives rise to the following inductive algorithm: 

Algorithm 4.3. The following algorithm computes the piecewise polynomial density of the 
push-forward of Lebesgue measure on M"q along {tj.) i— Yl^=i^k'^k- 

(1) Start with the unbounded chamber, where f = 0. 

(2) Iteratively jump over walls W separating the current chamber with an adjacent chamber: 



^In fact, the p-th summand of Theorem 4.1 is precisely the Duistermaat-Heckman measure corresponding to 
the isotropy representation of T on the symplectic vector space TpM, which is of course a non-compact symplectic 
manifold and, strictly speaking, does not fit into our setup. The decomposition of t* into regular chambers for th e 
moment map of M is refined by the common refinement of the chamber decompositions for the TpM (cf. 



il.ll 



^This is our reason for choosing a different definition for walls than the one used in jl2j. There, walls were 
defined as linear hyperplanes spanned by r — 1 linearly independent vectors. 
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(a) Denote by uji, . . . ,ujm the weights which lie on the h yperplane th rough W. 



(b) // the wall is minimal (m = r — 1), compute via Lemma 



c) Otherwise, recursively apply Algorithm 4-3 to compute the piecewise polynomial 







density of the push-forward of Lebesgue measure on M>q along (uk) i— X^fcLi ''J-k'^k 
On W itself it is given by a single homogeneous polynomial. Choose any polynomial 
extension to all of i* . 



(d) Compute the density on the adjacent chamber using (4.3). 



By combining Algorithm 4.3 with the Heckman formula, we arrive at the following algorithm 
for computing Duistermaat-Heckman measures. We shall call it the (Abelian) Heckman algorithm. 

Algorithm 4.4. Under the assumptions and using the notation of \Theorem 4-1 the following 



algorithm computes the piecewise polynomial density function of the Abelian Duistermaat-Heck- 
man measure: 

(1) Compute the density of each of the \M'^\ 
using \Algorithm 4-3 , 



iterated convolutions * Hcj^ ^ ★ . . . * He 



(2) Form their alternating sum according to Theorem 4-1 



The non- Abelian Duistermaat-Heckman measure can then be computed via Corollary 3.4- By 
passing to its support, we can also determine the non- Abelian moment polytope (cf. Remark 3.1). 



The algorithm as we have stated it assumes that the fixed-point data is part of the input. Let 
us describe it in the situations we are interested in: 

Remark 4.5. Consider the projective space M = P(y) associated with an arbitrary finite- 
dimensional, unitary /C-representation V . Torus fixed points in M correspond to weight vectors 
in V . Therefore, M"^ is finite if and only if all the weight spaces of V are one-dimensional. If this is 
the case, let V — ©^_q Cf^ be the weight-space decomposition, with Vf^ weight vectors of pairwise 
distinct weight cj^, so that the torus fixed points are precisely the points [vq], . . . , [vn] G M. Then, 
before renormalization, the isotropy weights in T[^^]M are given by the vectors ui — uj^ for / ^ k. 

Note that the representations associated with the pure-state quantum marginal problems 
displayed in [Table ~ indeed have one-dimensional weight spaces, so that Algorithm 4.4| is directly 
applicable: This is obvious for C^^ (g) . . . (8) and can also be verified for Sym (C^) and 
A^(C'^) (e.g., by observing that any single- row or single-column semistandard tableaux is already 
determined by its weight vector). However, other irreducible representations of SU(d), which 
correspond to indistinguishable particles of more exotic statistics, typically have weight spaces of 
dimension larger than one |33) . 

Remark 4.6. Consider more generally the action of T on a coadjoint i^-orbit M = induced 
by a group homomorphism ip: T ^ T ^ K. Even though this action might have infinitely many 
fixed points, there is an obvious way to write down an alternating sum formula for DH^, : Note 



that it follows directly from (1.2) that 



DH 



where tt = (d(/?)* is the dual map t* — )■ t*. Therefore, we can simply take the Abelian Heckman 
formula for the T-action (which is always applicable since the fixed point set of T is the Weyl 
orbit of A, hence finite), and push forward each summand along vr. In the case of a maximal- 
dimensional coadjoint orbit and for a suitable choice of renormalization direction, the result is 



just the push-forward of the Harish-Chandra formula (|3.1|) 



(4.4) 



with Qi, . . . , dj^ the positive roots of K. The formula for lower-dimensional coadjoint orbits can 
be obtained by using |10, Theorem 7.24] instead of (3.1). 



This density of course only depends on the hyperplane through W, and can therefore re-used for all other 
walls that span the same hyperplane. 
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In particular, this approach aUows the computation of the Abehan Duistermaat-Heckman 



measure for arbitrary setups of the quantum marginal problem by an obvious variant of Algo- 
Irithm 4.4| 



While Algorithm 4.4 and the variant described in Remark 4.6 are useful for computing low- 
dimensional examples, any approach relying on the Heckman formula has the major problem that 
the number of summands in the Heckman formula is typically very large (e.g., it is exponential 
in the number of distinguishable particles or fermions). Moreover, even though the Boysal- 
Vergne algorithm computes the density of a single summand chamber-by-chamber, this is less 
straightforward for the alternating sum, where all summands have to be evaluated in parallel. In 



4.2 below we will therefore derive an algorithm which does not suffer from these problems. 



There is also a non-Abelian Heckman formula due to Guillemin and Prato |39j (which suffers 
from the same problems). It can be deduced directly from the Abelian one by applying the 
derivative principle: 

Theorem 4.7 ( |39| (2.15)]). Suppose that there are only finitely many torus fixed points p G M"^ 
and that in each tangent space TpM each positive root a > or its negative occurs as an isotropy 
weight. Denote by rip the number of isotropy weights in TpM that are multiplied by —1 during 
renormalization. For each positive root a > and in each TpM, remove either a or —a from the 
list of renormalized isotropy weights. Denote the remaining weights by Wp^i, . . . ,ujp^n-R, o-nd let 
kp be the number of negative roots that have been removed. Then, 

DHf,= ^ (-l)"-+'=-5$^(p)^i/^^^*...*//<,^„_,|^^. 

In particular, the second assumption is satisfied when the moment map sends each torus fixed 
points to the interior of a Weyl chamber. 



Proof. Since dcjH±cj = i^o (cf- the proof of Proposition 3.1), the asserted formula follows at 



once by combining Corollary 3.4| with [Theorem 4.1 



Only the final remark needs elaboration: As observed by Guillemin and Prato, the assumption 
that ^k{p) S W ■ implies that the ET-stabilizer at each fixed point p is precisely T, so that 
the infinitesimal action of K generates a copy of t/t inside the tangent space TpM. Therefore, at 
any fixed point p, each positive root a > or its negative occurs as an isotropy weight. □ 



This gives rise to an obvious non-Abelian variant of Algorithm 4.4 



Algorithm 4.8. Under the assumptions and using the notation of \Theorem ^.7 the following 
algorithm computes the non-Abelian Duistermaat-Heckman measure: 

(1) Compute the \M'^ 



(see Remark 4-9). 



iterated convolutions ^^^ip ^ * . . .-kHcj^ using 



Algorithm 4-3 



(2) Form their alternating sum according to Theorem ^.7 



By passing to its support, we can also determine the non-Abelian moment polytope (cf. Re 



mark 3.1). 



Remark 4.9. There is a slight subtlety involved with the formulation of step (1) of Algorithm 4.8 
In case the renormalized isotropy weights Wp^i, . . . ,(2>p^n-B. in some TpM do not span all of t*, 
the corresponding iterated convolution is of course not absolutely continuous with respect to 



dX, and Algorithm 4.3 cannot be applied directly (see, e.g., the first proof of Proposition 5.1 ). 
Instead, we need to replace t* by the span of the a)p and apply Algorithm 4.3 accordingly. 



In g5 we will use both the Abelian and the non-Abelian version of the Heckman algorithm to 



compute the eigenvalue distribution of the reduced density matrices of a random pure state of 



two qubits (Proposition 5.1 



mixed states of two qubits (Proposition 5.8). 



and of bosonic qubits (Proposition 5.11 ), as well as of random 
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4.2. Single-Summand Algorithm for Projective Space. We will now derive explicit for- 
mulas for the Duistermaat-Heckman measure associated with a projective space, M = ¥{V), 
where V is a {n + l)-dimensional unitary representation of K, and where M is equipped with 
the Fubini-Study symplectic form wpSi normalized in such a way that its Liouville measure is 
equal to ^. The ET-action is Hamiltonian, and a canonical moment map is given by |58) 

(4.5) $^:P(F)^r, H^fx^-^^^ 

We start by decomposing the representation V into one- dimensional weight spaces, V = 
=0 Cufc, where Vk is a weight vector of weight Uk (repetitions allowed). In the corresponding 
homogeneous coordinates, the Abelian moment map has the following simple form, 

I 1 2 

(4.6) ^T:m)^i*, [^o:...:^n]^ ^f 'i",^' ^ 

z^fc=o Ffcl 

and it is straightforward to see that the Abelian Duistermaat-Heckman measure can be written 
as the push-forward of Lebesgue measure on the standard simplex along a linear map: 

Proposition 4.10. We have 

Here, P is the linear map M"+-'^ — t- t*, {tk) ^ X^fc ^k^k, o,nd dp is Lebesgue measure on the affine 
hyperplane H := {{tk) '■ "^Zk^k = 1} ^ M"^^, normalized in such a way that the standard simplex 
An := {{Pk) ■Pk>0, Ylk=oPk = 1} measure ^. 

Proof. The Fubini-Study measure is the push-forward of the usual round measure on the unit 
sphere ^^n+i ^ {(zq, . . . , Zn) ■ l^oP -|- . . . -|- \zn\'^ = 1} '^V along the quotient map {zq, . . . , z„) i— )• 
[zq '■ • • • : Zn], normalized to total volume ^. On the other hand, the round measure on the 
unit sphere also induces Lebesgue measure on the standard simplex by pushing forward along 
the map (zq, . . . , Zn) i— >• (l^oP, ■ ■ ■ , knP), since dxdy = d{r'^)d9 in polar coordinates. The claim 



follows from comparing (4.6) with P: {t^) i— )■ Yl^=Q^k^k- D 



Remark 4.11. Proposition 4.10 can also be established by applying the Heckman formula as 
described in [Remark 4.61 

Denote by dp/d\ a differential form corresponding to Lebesgue measure on the affine subspaces 
P^'^{\) n H, normalized in such a way that 

(4.7) dp = dp/dXAP*{dX) 

when restricted to the affine hyperplane H. 

Proposition 4.12. The density function / it*—)- [0, oo) of the Abelian Duistermaat-Heckman 
measure is given by 

n n 



/(A) = vol{pk > : ^Pfe^A, = A,^pfc = 1}, 



k=0 k=0 



where the volume is measured with respect to the measure induced by dp/dX on P (A) Pi H. 
Proof. For all test functions g G Cb{i*), we have 



(DH^(^),5)=/ dpg{P{p))= f dxi f dp/dx]g{X), 
Ja„ Jt* Wp-i(A)nA„ / 



by using (4.7) and Fubini's theorem for the fibration |41[ pp. 307]. □ 

That is, the Abelian Duistermaat-Heckman density measures the volume of a family of convex 
polytopes parametrized by t*. This is also true for the density of the iterated convolutions 



studied in §4.1| (see (4.9) below). There are exact numerical schemes that can be used to 



compute the polynomial density functions on each regular chamber which have already been 
implemented in software packages, e.g., the parametric extension of Barvinok's algorithm t4j 



described in [lOH 1100] . We will not pursue this route any further. However, in §6.3 we will show 
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that its "quantized" counterpart gives rise to an efficient way of computing the corresponding 
representation-theoretic quantities (in particular, the Kronecker coefficients). 

In the foUowing, we wiU instead describe a combinatorial algorithm based on the same principles 
as our Heckman algorithm. Before doing so, let us determine explicitly the regular chambers 
for the Abelian moment map, i.e., the connected components of the set of regular values of $t, 
each on which the measure is given by a polynomial. For this, we define the support of a point 
p = [v] £ ¥{V) as the set of weights which contribute to the weight-space decomposition of v, 

suppp := {ujk : Zk ^ 0,p = [zo : . . . : Zn]}- 

The significance of this definition is that the support of a point already fully determines whether 
it is regular or singular: 

Lemma 4.13. Let p G P{y)- Then p is a regular point of the Abelian moment map if and only 

spanjo; — u : uj , uj' G suppp} = t*. 

Proof. It follows readily from the definition of the moment map that a point p is regular if and 
only if tp, the Lie algebra of its stabilizer, is trivial {42| Lemma 2.1]. But ip is already determined 
by the support of p: 

tp = {X €t: uj{X) = uj'{X) Vw, Lo' e suppp} 
This is the annihilator of the linear span in the statement of the lemma. □ 

We arrive at the following characterization of the set of singular values of the Abelian moment 
map: 

Proposition 4.14. The set of singular values of is the union of all convex hulls of subsets 
containing (at most) r weights, 

convjwfc : k £ 1} = convjwfc : G /}. 

#I=r #I<r 



Proof. It is clear from (4.6) and Lemma 4.13 that the convex hull of any subset of weights of 
cardinality at most r consists of singular values. The converse follows from Caratheodory's 
theorem. □ 

In particular, the singular walls are convex hulls of r weights in general position. From this 
description we can easily determine the regular chambers. Observe again that there is a single 
unbounded regular chamber. 



We will now use the result of Boysal and Vergne described in §4.1 to derive intrinsic formulas 



for the jumps of the Duistermaat-Heckman density when crossing a singular wall. Recall that 
the measures they consider are push- forwards of Lebesgue measure on the convex cone M>q^ 
rather than of Lebesgue measure on the standard simplex A„, which is of course the intersection 
of M.y'^^ with the affine hyperplane H = {(t/t) : Yl^=o^k = !}■ It is however straightforward to 
translate between both pictures: In order to avoid confusion, we shall use the same convention 
as in §4.1 that hatted quantities correspond to the Boysal- Vergne picture. Let us consider the 
"extended" weights cD^ := (u^, 1) St*©M {k = 0, . . . ,n) together with the corresponding linear 
map 

n n 

P: M"+i ^t*©M, {tk)^J2^kC0k = {P{to,...,tn),^tk). 

k=0 k=0 
Denote by dt standard Lebesgue measure on M""*"^ and equip t* ® M with the measure dX = dXds, 
where ds is standard Lebesgue measure on M. Choose a differential form dt/dX inducing Lebesgue 
measure on the fibers of P, normalized in such a way that 

(4.8) dt = dt/dX A P*{dX) = dt/dX A P*{dX) A {dto + ... + dt^). 
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Then one can establish just as in the proof of [Proposition 4.12 the fohowing formula for the 
density function of the push-forward of Lebesgue measure on M"g^ along P with respect to 

dX = dXds, 

n n 

(4.9) /(A, s) = vol {tfe > : J] t^uju = A, J]] tfc = s}, 



fc=0 



fc=0 



where the volume is measured with respect to dt/dX. But comparing (4.7) and (4.8) and noting 
that dt = dp f\ {dto + . . . -\- dt^) on H, we see that in fact dt/dX and dp/dX induce the same 
measure on the fibers P~^{X) = P^^{X, 1), so that 



(4.10) 



DH 



;^(^) = /(A) dX = fix, 1) dX. 



This shows that we can work equivalently in the convex cone picture of Boysal and Vergnej^ 

We shall now describe the jump formula. Let Vl^ be a singular wall separating regular chambers 
A± C t*, and choose a normal vector ^ S t* pointing from A_ to A_|_. Order the weights such that 
precisely uq, . . . ,uJm-i lie on W. Denote by dw Lebesgue measure on the hyperplanes parallel to 
W, normalized in such a way that 

(4.11) dX = dwAd(„ 

where d^ is the pullback of the standard volume form of M along the coordinate function { — ,£,)■ 
Denote by f± the polynomials describing the density function / on the regular chambers A±. 
Finally, consider the Duistermaat-Heckman measure for the action of T on the projective space 
over V\Y = ®^!lo^ Cffc, the direct sum of the weight spaces corresponding to the weights which 
lie on the hyperplane through W. Its density with respect to dw is given by a single polynomial 
on the singular wall W, since W is always contained in the closure of a regular chamber for 
IF'(Hv)- Choose any polynomial function fw extending it to all of t*. 

Proposition 4.15. The jump of the Abelian Duistermaat-Heckman density across the singular 
wall is given by 

f+W - /-(A) = Res I /iy(9x,9y)™ TT—^ ^— 

iz-0 y w Ylj^^^ z{uk - wo, Cl + K, x) + y j ^^^^^^^ 

Here, f-^{X,s) = s'^~'^ fw{^) is the homogeneous "extension" of fw to i* © M. 



Proof. The convex cones A± through A± x {1} are chambers in the sense of Boysal and 
Vergne. They are separated by a wall W, namely the convex cone through W x {1}. Note 
that C = (C) ~ (^OjO) is ^ normal vector to W. Denote by f± the homogeneous polynomials 
describing the density function of the push-forward of Lebesgue measure on M"g^ along P. It is 

clear that dw = ds A dw induces Lebesgue measure on W and that it is normalized in such a 
way that dX = dw A d(,. By (4.10) and the jump formula (4.3) of Boysal and Vergne, we have 

/ ^ z((X,l),x+zi) \ 

/+(A) - /-(A) = /+(A,1) - /_(A,1) = Resl^^g /^(5,)— — — ^ 

The polynomial f^ as defined above agrees with its original definition in §4.1 since it is a 

[y , which b 
(x, y) G t* 



homogeneous polynomial and can thus be reconstructed from fiy, which by (4.10) is its restriction 
to the slice t* x {1}, by the formula given above. Writing x = {x,y) G t* M and expanding the 
hatted quantities, we arrive at the assertion. □ 



As in §4.1 the case where only a minimal number of weights lie on the affine hyperplane 
through W is particularly simple to evaluate: 

'^The push-forward of Lebesgue measure on R"^'^ along P can also be understood as the Duistermaat-Heckman 
measure associated with the Hamilt onian T x U(l)-action on the complex vector space V, where U(l) acts by 



scalar multiplication (cf. Footnote 1 1 
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Lemma 4.16. Suppose that precisely r weights ojq, . . . ,ujr-i He on the affine hyperplane through 
W. Then, 



fw^ = f\v = \dX [ uji - u;o, . . . , ujr-i - Wo, 



Proof. We argue as in the proof of Lemma 4.2 In view of Proposition 4.10 and the minimahty 
assumption, the map (qk) i— >• Yl^k=o Ik^k along which we push forward is an isomorphism, and 
fw is equal to the constant of proportionahty between the push-forward of Lebesgue measure on 
H (normalized in such a way that the standard simplex has measure ^) and the measure dw. 
We can compute this constant by comparing the volume of the parallelotope spanned by the 
(uJk)' For the former measure this constant is one (by its very normalization), while for the latter 
it follows from ( 4.11[ ) that 

lief. 



dX{u}i — UJQ, . . . , UJr-l — OOq, ^) = duj{uJi — bjQ, . . . ,ljJr 



□ 



These results give rise to the following inductive algorithm for computing the Abelian and 
non-Abelian Duistermaat-Heckman measure of a projective space. We will call it the single- 



summand algorithm, since in view of (4.10) it amounts to computing a push- forward measure 



that is equivalent to a single summand of the Abelian Heckman formula (cf. Theorem 4.1). 



Algorithm 4.17. The following algorithm computes the piecewise polynomial density function 
of the Abelian Duistermaat-Heckman measure of the projective space ¥{V): 

(1) Start with the unbounded regular chamber, where / = 0. 

(2) Iteratively jump over singular walls W separating the current regular chamber with an 
adjacent regular chamber: 

(a) Denote by CjJq, . . . , w^-i the weights which lie on the hyperplane through W . 

(b) // the wall is minimal (m = r), compute fw via Lemma 4-16 



(c) Otherwise, recursively apply Algorithm ^.i? to compute the piecewise polynomial den- 
sity of the Abelian Duistermaat-Heckman measure ofWiVw), where Vw = ®T=o '^'^k 
is the direct sum of the weight spaces for the weights in (a)^ On W itself, it is given 
by a single polynomial. Choose any polynomial extension fw to all of t* . 



(d) Compute the density on the adjacent chamber using Proposition 4-15 



The non-Abelian Duistermaat-Heckman measure can then be computed via \Corollary 3.4 ' By 
passing to its support, we can also determine the non-Abelian moment polytope (cf. Remark 3.1). 



Remark 4.18. In view of (2.5) and (2.6) and by passing to the purified double (§2.3), Algo- 



rithm 4. 171 solves the problem of computing the eigenvalue distribution of reduced density matrices 
in complete generality. 

We conclude this section by explicitly stating the Abelian and non-Abelian jump formula for 
the case where only a minimal number of weights lie on the affine hyperplane through the wall. 
They will be used later for computing examples. 

Corollary 4.19. Suppose that precisely r weights wq, . . . , (^r-i lie on the affine hyperplane 
through the singular wall W . Then the jump of the Abelian Duistermaat-Heckman density across 
the wall is given by 

' (A-a;o,0"~' 



f+{X)-fA\) = fw[\{{i^k-^o,i) 



\k=r 



(n — r)\ 



where fw is the constant from Lemma ^.ifi} 

Proof. This follows immediately from [Proposition 4.15 by pulling out the constant fw, setting 
X = y = and evaluating the residue at z = 0. □ 



The non-Abelian formula follows directly by applying [Corollary 3.4 



This density of course only depends on the hyperplane through W, and can therefore be re-used for all other 
singular walls that lie on the same hyperplane. 
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Figure 2. Abelian moment polytope of two qubits (grey square) and its de- 
composition into four bounded regular chambers by the singular walls (thick 
lines). 



Corollary 4.20. Suppose that precisely r weights ojq, . . . ,ujr-i lie on the affine hyperplane 
through the singular wall W , and that n — r > R, so that the non-Abelian Duistermaat-Heckman 
measure ofF(y) is absolutely continuous in the vicinity ofW. Denote by the polynomials 
describing its density on the regular chambers. Then the jump across the wall is given by 



/f (A) - /i^(A) = fw n^^fc -^O'O n 



in 



Ry. 



where fw is the constant from Lemma 4-16 



Remark 4.21. Corollary 4.19 has already been established in |37j . where the authors also envisaged 



an algorithm similar to our Heckman algorithm. They did however not have a general jump 



formula such as (4.3) at their avail. Instead, they had to resort to an inexact formula which in 



general only holds in highest order (in the distance to the wall). 



5. Examples 

In this section we illustrate our algorithms by computing some eigenvalue distributions of 
reduced density matrices. The global quantum states will always be chosen according to one of 
the invariant probability measures described in §2.1 Many of our examples will involve qubits, 



i.e., quantum systems modeled by two-dimensional Hilbert spaces, so that the algorithms can 
be nicely visualized. But of course our algorithms can be used to determine the eigenvalue 



distributions for arbitrary instances of the quantum marginal problem (see Remark 4.18). 



5.1. Pure States of Multiple Qubits. We start by considerin g pur e states of N qubits, where 



by tensor products (cf. §2.2). It will be convenient to 



K = SU(2)^ acts on M 

identify t* = in such a way that the positive Weyl chamber corresponds to the cone M>q 
and the fundamental weights to the standard basis vectors ej = {5j^k) (^ = 1, • • • That 
is, if A = (A^-'-') € t* then we will by slight abuse of notation identify A'--'^ with the scalar 



i{X^P — ^2^) = 2iX^f' . It follows that dX is simply the usual Lebesgue measure on M^, that the 
symplectic volu me po lynomial is given by Pk{X) = A^^^ • • • X^^\ and that the positive roots are 



2ei, . . . , 2e]sr (cf. §1.1 and (2.3)). Moreover, (2.2) amounts to assigning to a point (A^^)) G the 



tuple {pi, . . . , pn) of diagonal density matrices acting on C , where pj has maximal eigenvalue 



A 



l+iX[ 



U) _ i+A(i) 



max — 2 ^ — 2 

We first discuss in detail the toy example of A'^ = 2 qubits, demonstrating both the non-Abelian 
Heckman algorithm and the single-summand algorithm. 

Proposition 5.1. The non-Abelian Duistermaat- Heckman measure for the action o/SU(2) x 
SU(2) on ® £?) is given by 

jSU(2)xSU(2) », 1 



(DH 



P(C2®C2) 



f{t,t)dt, 



EIGENVALUE DISTRIBUTIONS OF REDUCED DENSITY MATRICES 



25 



i.e., by a one- dimensional Lebesgue measure supported on the diagonal between the origin and 
(1,1)- 



Proof using the non-Abelian Heckman algorithm {Algorithm 4-S). The four fixed points of the 
action correspond to the standard basis vectors Cj <^ (j, k = 1,2), which are weight vectors 
of weight (ibl,ibl) using the conventions fixed above (the vertices of the grey rectangle in 
Figure 2^ . Let us choose the direction 7 = (—2,-1) for renormalization. After removal of 
the positive and negative roots, (±2,0) and (0, ±2), only a single renormalized isotropy weight 
remains at each fixed point (cf. [Remark 4.5 ). Therefore, Theorem 4.7 shows that the non-Abelian 
Duistermaat-Heckman measure is given by the restriction to the positive Weyl chamber of 

^{1,1} * -f^(-2,-2) - '^(1,-1) * H{-2,2) + ^(-1,1) * H{-2,2) - '^(-1,-1) * -f^(-2,-2)- 

Only the first summand contributes to the positive Weyl chamber, and its restriction is given 
precisely by the formula displayed above (cf. Remark 4.9). □ 



Proof using the single- summand algorithm (Algorithm Note that the non-Abelian wall 

jump formula (Corollary 4.20) is not directly applicable, since n — r ^ R. Indeed, as we have 



seen above, the non-Abelian measure does not have a Lebesgue density, since it is concentrated 
on the diagonal. 



Therefore, we will follow Algorithm 4.17, which uses the Abelian wall jump formula, and 
afterwards takes partial derivatives in direction of the negative roots according to [Corollary 3.4 



The decomposition of t* into regular chambers is indicated in Figure 2 We start in the unbounded 
chamber, where the density is equal to the zero polynomial and cross the horizontal singular 
wall at the top. Evaluating the Abelian jump formula (Corollary 4.19 say, with lvq = (1, 1) 
and ^ = (0,-1)), we find that the density polynomial on the upper regular chamber is equal to 

Next, we cross the diagonal singular wall separating the upper and the right-hand side regular 
chamber. Using the Abelian jump formula once again, we see that the density polynomial changes 

by 1(A(2)-A«). 

Therefore, the Abelian Duistermaat-Heckman measure has the following piecewise polynomial 
density on the positive Weyl chamber: 



1 



max 



(A«,A(2)) 



Taking partial derivatives in the direction of the negative roots, (—2, 0) and (0, —2), we arrive at 
the measure asserted above. □ 

Corollary 5.2. The joint distribution Pdg of the maximal eigenvalues of the reduced density 
matrices of a randomly-chosen pure quantum state of two qubits is given by 



(Peig,/) = 24 



f{s,s) ( s - ^ ) ds, 



for all test functions /(AmL, AmL) 



Proof. According to (2.5), multiply the non-Abelian Duistermaat-Heckman measure by the 
symplectic volume polynomial Psu{2)xSU(2)(A) = X^^^X^'^\ divide by ^, the volume of P(C^ C^C^). 

Finally, push forward along (A^^)) ^ (A^i?L = ^^^). □ 

This eigenvalue distribution is in fact known more generally for bipartite pure states chosen at 
random |76l 1105) . We will later show how to compute its generalization using the techniques of 
this paper (Corollary 5.15). 



For higher tensor powers, evaluating the Heckman formula quickly becomes unwieldy. However, 
it can still be used to compute the Duistermaat-Heckman measure locally: 
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Figure 3. Non-Abelian moment polytope of three qubits and its decomposition 
into six bounded regular chambers. 



Proposition 5.3. The non-Abelian Duistermaat-Heckman measure for the action o/SU(2)^ on 
is on the closures of the regular chambers that contain the vertex (1, . . . , 1) given by 
the convolution product 

where {uJk} is the set of weights of the form (—2, . . . , —2, 0, . . . , 0) (at least two non-zero entries) 
as well as their -permutations. 

Proof. If we renormaUze with respect to the direction 7 w (—1,...,—!) then just as in the 



first proof of Proposition 5.1 only a single summand in the non-Abelian Heckman formula 
contributes in the vicinity of the vertex (1, . . . , 1) and, moreover, this summand is of the above 
for m: Indeed, th e weights {ook} are precisely the isotropy weights with the negative roots removed 
(cf. Remark 4.5). Since the density function of DH^ is polynomial on each regular chamber 



adjacent to the vertex, we can extend the local formula to their closures. □ 

It is in fact easy to see that the domain of validity of this formula is the intersection of the 
half-space 

N 

A : ^^^^ >N-2 
i=i 

with the positive Weyl chamber (the regular chambers not adjacent to (1, . . . , 1) lie in the 
complement of this half-space) . 

Remark 5.4. [Proposition 5.3| gives a local description of the non-Abelian moment polytope, 
namely by the cone based at (1, . . . , 1) and spanned by the rays with direction vectors {uJk}- 
By convexity, its intersection with the positive Weyl chamber is an outer approximation to the 
moment polytope. 

Let us specialize to the case = 3: Here, precisely the rays with the direction vectors 
(—2, —2, 0), (—2, 0, —2) and (0, —2, —2) are extremal. Their intersection with the positive Weyl 
chamber has to be contained in the non-Abelian moment polytope: Otherwise, there would 
be additional vertices in the interior of the positive Weyl chamber — but only (1, 1, 1) is the 
image of a torus fixed point. Since also the origin is contained in the moment polytope (the 
Greenberger-Horne-Zeilinger state, [ip] = [ei <^ ei ei + €2 62], is a preimage of the origin 
[36j). we conclude that the convex hull 



conv{(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 1)} 
is an inner approximation to the moment polytope. Both approximations are in fact equal and 



therefore describe the moment polytope precisely (Figure 3). Inequalities characterizing the 



moment polytope for A^ qubits have been determined in [53j. 
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Proposition 5.5. The non-Abelian Duistermaat-Heckman measure for the action o/ SU(2)^ on 
$5 (X" C^) has the piecewise linear Lebesgue density 



( J_ 

16 

J_ 
32 



min A^-'^ 



in the lower pyramid, 
mm A^-') ] %n the upper pyramid, 
otherwise 



(compare Figure 3). 



Proof. By [Proposition !x3 the non-Abelian Duistermaat-Heckman measure is on the closures of 
the regular chambers containing (1,1,1) given by the convolution 



'^(1,1,1) * H{-2,~2-2) * -f^{-2,-2,0) * ^(-2,0,-2) * -^'(0,-2, -2) • 



Using (4.1) we can readily compute its density: 



dti 



dU6{(l)+t,[ll)+t,[--l^+t,[-l)+t,[-'2 



A) 



32 



f°° ( 1 

J — oo ^ 

l^"d,,le((l-.0(l)-A) 

1 ^ 

— max{0,l -^A(^) + 2min(A(^))}, 



A) 



where Ic is the indicator function of the cone spanned by ( i 



0^, and (^1 i-e., 



1 if a + 6 > c and a + c > 6 and b + c> a, 
otherwise. 



lc{a,b, c) 

We have therefore established the claimed density on the complement of the lower pyramid. 



According to Corollary 4.20, the jump across the hyperplane separating the upper and the 



lower pyramid is given by 



K-a-('-i")=^(t 



A(^-) - 1 



(the left-hand side terms are ordered just like in the jump formula). This is precisely the 
difference between the densities on the upper and lower pyramid as asserted in the statement of 
the proposition. □ 



It is straightforward to deduce from this the eigenvalue distribution (cf. the proof of Corol- 



lary 5.2): 



Corollary 5.6. The joint distribution of the maximal eigenvalues of the reduced density matrices 
of a randomly-chosen pure quantum state of three qubits has Lebesgue density 



8! n ^ 




min Amax - o 



max 



{o4(i-E?.iA 



^max I ~l" min Amax 



if A'^-'^ < 2 

y Z^j=l ^max ^ ^1 

} ^/E-=iA^j^Ix>2, 



on the space of maximal eigenvalues (Amix) G [\, 1]^- 



Remark 5.7. Our use of the local convolution formula (Proposition 5.3) and of the non-Abelian 



wall jump formula (Corollary 4.20) were merely convenient shortcuts: It is clear that we could 



have completely algorithmically computed the measure by following Algorithm 4.17 
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Figure 4. (a) Chambers and support (gray) and (b) density function of the 
iterated convolution computed in Lemma 5.9[ 



5.2. Mixed States of Two Qubits. We wih now use the non-Abehan Heckman algorithm to 
treat the case of random two-qubit states with fixed, non-degenerate global eigenvalue spectrum. 
That is, we consider the action oi K = SU(2) x SU(2) on a coadjoint SU(4)-orbit through a point 
A contained in the interior of the positive Weyl chamber. 

Recall that the Weyl group of SU(4) is the symmetric group S4, with equal to the 



signum of a permutation w G S4. By (4.4), 
(5.1) 



DH 



where ai, . . . ,aQ are the positive roots of SU(4) (see §1.1 for our conventions), and where vr is 
the restriction map t* — t- t*, with t* the dual of th e Lie algebra of the maximal torus of SU(4). 
With respect to our identification t* = M? fixed in §5.1 , the map vr is given by 

(5.2) 7r:t*^M2^ (Ai, . . . , A4) ^ 2i(Ai + A2, Ai + A3). 

One computes readily that the — 7r(afc) are precisely the weights (—2, 2), (—2, 0) (twice), (—2, —2) 
and (0, —2) (twice). In particular, the two negative roots of SU(2) x SU(2) are contained in this 
list (each of them is in fact contained twice) . By applying [Corollary 3.4 we arrive at the following 
formula: 

Proposition 5.8. The non-Ahelian Duistermaat-Heckman measure for the action o/SU(2) x 
SU(2) on a coadjoint SU(4) -orfcii Oy^ with A € tJ^Q is given by 



DH 



SU(2)xSU{2) 



* H{-2,2) * ^(-2,0) * -f^(-2,-2) * ^(0,-2) 



Following Algorithm 4.8 we evaluate the right-hand side iterated convolution using Algo- 



rithm 4.3 The result is the following: 

Lemma 5.9. The measure ^^(_2,2) *-f^(-2,o) * -^^(-2,-2) *-f^{o,-2) ^^-^ Lebesgue density 

in chamber 0, 



/(A«,A(2)) 



^ (A(^) + A*^^)) in chamber 1, 

^ ((A(i))' + 2A(i)A(2) - (a(2))2) tn chamber 2, 



in chamber 3. 

See Figure ^ for the labelling of the chambers and an illustration of the density. 

The density of the non-Abelian Duistermaat-Heckman measure is thus given by the restriction 
to the posit ive Weyl chamber of an alternating sum of 24 copies of the density described in 
Lemma 5.9 one copy attached to each of the points 7r(?I;A). In view of the geometry of the 
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support of the latter density, it is clear that in fact only summands for points in the right 
halfplane {A^^^ > 0} contribute (i.e., at most half of the points). Using (5.2), one finds that the 



points tt{w\) are the six points whose coordinates are equal to any two out of the three values 
ci = |2i(Ai + A4)|, C2 = 2z(Ai + A3), or C3 = 2z(Ai + A2) (without repetitions), as well as their 



Weyl conjugates. See Figure 5] for illustration. 



Moreover, we can deduce that the non-Abelian moment polytope has the form described in 
the figure. To do so, we simply need to check in each regular chamber whether the density 
polynomial vanishes. By doing so and describing the resulting polytope in terms of inequalities, 
we recover a well-known result by Bravyi [13j: 

Corollary 5.10. The non-Abelian moment polytope for the action o/SU(2) x SU(2) on a generic 
coadjoint S\J (4:) -orbit is given by 

Asu{2)xSU(2)(Oa) = {(A('\a(2)) : < A«,A(2) < c3,A« + A(2) < ^2 + C3, |A« - A^^)] < C3-C1}. 
In the limit where the global state becomes pure, the moment polytope converges to the 



diagonal between the origin and (1, 1). This is in agreement with Proposition 5.1 One can 
similarly recover the eigenvalue distribution of the reduced density matrices of a random pure 
state of two qubits by taking a corresponding limit. 



In view of Proposition 2.8 the distributions computed in Proposition 5.8 can be assembled to 
give the joint eigenvalue distribution of the reduced density matrices of a randomly-chosen pure 
state in F(C^ ^ ^ C^). 



5.3. Pure States of Bosonic Qubits. We now turn to random pure states of bosonic qubits, 
where K = SU(2) and M = P(Sym^(C2)). We wiU use the Abelian Heckman algorithm: 

Proposition 5.11. The Abelian Duistermaat-Heckman measure /or P(Sym^(C2)) has Lebesgue 
density 



^ ' k=-N,~N+2,...,N 



_ ^r_l_^ m \ 2 / 



N- 



Here, we set (A — k)^~^ = (A — k)^~^ for \ > k and otherwise. 

Equivalently, N\ times the Abelian Duistermaat-Heckman measure is equal to the probability 
distribution of the sum of N independent random variables that are uniformly distributed on the 
interval [—1,1] [31, §1.9, Theorem la]. 

Proof. The weights of Sym^(C'^) are {—N, —N -\- 2, . . . , N}; let us write Vk for a weight vector 
of weight k. The associated projective space has precisely N -\- 1 torus fixed points. At any such 
fixed point [vk], the isotropy weights are given by 

{{l-k) -.1 = -N,-N + 2,...,k,...,N}, 
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Figure 6. (a) and (b) Density of the Abelian and the non-Abehan Duistermaat- 
Heckman measure for P(Sym^(C^)), and (c) corresponding maximal eigenvalue 
distribution of the one-body reduced density matrix. 



and we will denote them by cj^ i^, . . . ,Lbk,N (cf- Remark 4.5). Observe that precisely = ^^^^ 
of them are negative with respect to the renormalization direction 7 = +1. By (4.1), the 
corresponding summand of the Heckman formula is equal to the push-forward of Lebesgue 
measure on M^q along 

N 

P: M>o u*, (si, ...,sn)^ '^Si\Lbk,i\ + k. 
We first compute its cumulative distribution function: 



i=l 



P^{ds) ((-00, {k + A)]) = ds({si, . . . , SAT > : V'sjlwfc.il < A}) 



TV 

1=1 

N 



ds{{si, 
1 



, SN 



>0: J^s, <!}) 



m 2^(^)!(^)! 
1 / N 



i=l 

Ai^ 



N 



The density is then given by the derivative, 

1 



and by forming the alternating sum of these terms we arrive at the formula displayed above. □ 



N 



In §7.2 we give an alternative proof of Proposition 5.11 using representation theory and 
combinatorics. It follows directly from the derivative principle that the non-Abelian Duister- 
maat-Heckman measure is given by the following formula: 

Corollary 5.12. The non-Abelian Duistermaat- Heckman measure for the action o/SU(2) on 
P(Sym^(C^)) with N > 2 has Lebesgue density 



2^-i(A^-2)!A^! 



iV+fc 
2 



-N-N+2,...,N 



on [0, 00). 

Again, it is clear how to translate the above into the eigenvalue distribution of the one-body 
reduced density matrix by using (2.6). See Figure 6 for an illustration in the case of = 2 
bosonic qubits. 

As an application, let us compute the average value of the reduced purity of a randomly-chosen 
pure state of bosonic qubits. The reduced purity is by definition equal to 

v2 I I'l \ \2 



(5.3) 



IP1II2 



'^max + (1 - Ar 



where pi denotes the one-body reduced density matrix and Amax its maximal eigenvalue. 
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Corollary 5.13. The average reduced purity (5.3) of a randomly-chosen pure state of N bosonic 
qubits is given by 

1 1 

2 ^ 2N' 



Proof. We will not use Corollary 5.12 directly but instead work with the Abelian Duistermaat- 



Heckman measure: Denote by P the probability distribution of the maximal eigenvalue of the 
one-body reduced density matrix. By using (2.6) and Amax = 5 + 2^1 ^^'^ that the average 
reduced purity is given by 



A 



+ 1 - A„ 



dP(A„ 



1 



4Ar2 
1 



X (^{N + Xf + {N - Xf^ Nl 



dDR 



SU(2) 
M 



(A) 



(A). 



By the derivative principle, Corollary 3.4 this is equal to 



1 

iV2 
1 



(iV^ + 3A^) NldDR 



U(l) 
M 



(A) 



2iV2 



[N^ + 3X^) N\dBRlj^\x). 



In [Proposition 5.11 we have seen that A^! DHYr^^ is the probability distribution of the sum of N 



independent random variables that are uniformly distributed on the interval [—1,1]. Since the 
variance of any such random variable is | and since variances of independent random variables 
are additive, the above is equal to 



1 



2iV2 



N 



1 1 

2 ^ 2N' 



□ 



In accordance with the concentration of measure phenomenon, pi — )• 1/2 in distribution as 
— 7- 00. We remark that our result matches |8H Theorem 34] if one works out the quantities 
left uncalculated therein. Note that our proof illustrates the power of the derivative principle: 
Instead of explicitly computing the eigenvalue distribution, we can reduce to the Abelian Duister- 
maat-Heckman measure by differentiating the quantity we are interested in. 



5.4. Pure States of Bipartite Systems. We conclude this series of examples by re-deriving 
the eigenvalue distribution of the reduced density matrices of a randomly-chosen pure state in 
the case of a general bipartite quantum system, corresponding to the action of SU(a) x SU(6) 
on M = P(C C^). Instead of following one of the algorithms it will be most convenient to 
directly work with the formula given in [Proposition 4.12 



Suppose that b > a. If6>o-|-l then Assumption 1.1 is not satisfied (Lemma 2.9): Indeed, it 



always follows from the singular value decomposition that 



(5.4) 



eigPB = {eigpA,0, • • • ,0), 



so that in this case the non- Abelian moment polytope is contained in the boundary of the positive 



Weyl chamber. However, (5.4) of course implies that for any choice of 6 > a the joint eigenvalue 



distribution is already determined by the eigenvalue distribution of pA- That is, it suffices to 
compute the Duistermaat-Heckman measure for the action K = SU(a). Denote by T the 
standard maximal torus of SU(a). As in (2.1 ), we identify points A G t = t* with diagonal density 
matrices A = ^ -|- iA. 

Clearly, the Abelian moment polytope consists of those A with A G Aa_i, and the non- Abelian 



moment polytope is its intersection with the positive Weyl chamber (cf. §2.3). 
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Lemma 5.14. On the Ahelian moment polytope, the Duistermaat-Heckman measure for the 
action of the maximal torus of SU(a) is proportional to 

j=l j=i V / 

Proof. Choose the weight-space decomposition of (8) given by the standard basis vectors 
Cj ® Ck (j = 1, . . . , a and /c = 1, . . . , 6). According to Proposition 4.12 the density of the Abelian 
Duistermaat-Heckman measure is given by 

/(A) = vol{pi,i, . . . ,pa,h > : ^Pi,fc = Xj = - + iXj (j = 1, . . . ,o)} 

fc=i " 

with respect to the volume measure dp/dX defined therein. Note that the right-hand side set is 
the Cartesian product of a rescaled standard simplices. The measure factorizes accordingly, and 
it is easy to see that 

a b 1 

/(A) = l[vol{pi,...,Pb>0:^Pk = Xj} = -l[ Aj-i 
j=i k=l j=i 

for A G Aa_i, and zero otherwise, with Z a suitable normalization constant. □ 

Corollary 5.15. On the non- Abelian moment polytope, the Duistermaat-Heckman measure for 
the S\J {a)-action on F{C"' (g) C^) is proportional to 

n ^'''^ n cxj - Xk) dx. 

j=l j<k<a 

Proof According to the derivative principle, we have to apply nj<A: ^i'^^k ~^>^i ) ~ Wj<ki^\f, ) 
to the Abelian Duistermaat-Heckman density as computed in [Lemma 5.14| 

This is a partial differential operator of order (2) , therefore the resulting non- Abelian density 
polynomial has total degree at most dmax = ^(^ ~ 1) ~ ^(^ ~ Since we differentiate each 

variable at most a — 1 times, it is a multiple of the symmetric polynomial 11^=1 
other hand, the result is evidently antisymmetric, and therefore a multiple of the Vandermonde 
determinant nj<fc('^j ~Xk)- Since the total degrees add up to dmax) this implies the assertion. □ 



In view of (2.5), this result implies the following well-known formula |76 l 1105) : 

Corollary 5.16. The distribution of the eigenvalue spectrum X = eig/?i of a randomly- chosen 
bipartite pure state p on C"" C^) has Lebesgue density proportional to 

fm-^ n (^.-^^)^ 

j=l j<k<a 

on the space of eigenvalue spectra {X € Aa-i : Ai > . . . > Aq}. 

It is also easy to deduce the corresponding formula for the action of SU(a) x SU(6): 

Corollary 5.17. Denote 6?/ A = {A € : A G Aa-i} the non-Abelian moment polytope for 
the S\J {a)-action. Then the push-forward of Liouville measure along the moment map for the 
SU(a) X S\J{b)-action is given by 

>• {{A,0,...,0)-i)/i 



where Z is a suitable normalization constant. 
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Proof. By (5.4), each coadjoint orbit Ox for SU(a) is paired with the coadjoint SU(6)-orbit 
through ((A, 0, . . . , 0) - By its SU(a) x SU(6) -invariance, on each such pair of coadjoint 

orbits the push-forward measure is just a multiple of the usual Liouville measure. The assertion 



follows by observing that the density in Corollary 5.15|is at any point A proportional to the 



symplectic volume of the corresponding coadjoint SU(6)-orbit. □ 

6. Multiplicities of Representations 

All results discussed so far can be considered as asymptotic limits of corresponding statements 
in representation theory, at least if the Hamiltonian iC-manifold M can be linearized ("quantized") 
in a certain technical sense. This is in the spirit of Kirillov's orbit method and the theory of 
geometric quantization fUl |3H 11031 EH El] . 

In particular, this is the case when M is a if-invariant smooth irreducible complex projective 
subvariety of ¥{V) for a finite-dimensional unitary i^-representation V. In this situation, the 
Fubini-Study form of ¥{V) restricts to a non-degenerate symplectic form on M, and the ii'-action 



is Hamiltonian with moment map the restriction of (4.5). We still assume that Assumption 1.1 
is satisfied. 

Coadjoint orbits Ox through dominant integral weights A G A* (and only these) can be realized 
in this setup [56j: They are in a natural way projective subvarieties of ¥{Vx), where Vx is the 
unitary i^-representation with highest weight A. 

In particular, the quantum marginal problem can be analyzed in this framework: Coadjoint 
SU((i)-orbits through integral highest weights correspond to Hermitian matrices with integral 
eigenvalue spectra, and it suffices to consider these, since we can always rescale and take limits. 



or simply pass to the purified double (§2.3) 



In §6.1 we will recall the limit alluded to above. We then proceed to describe the representation- 



theoretic analogue of the derivative principle: Multiplicities of irreducible iC-representations can 



be computed from weight multiplicities by taking finite differences (Lemma 6.2). In the case of 



the projective space associated with a unitary ^-representation, the relevant weight multiplicities 



are those for the symmetric powers of the representation. In §6.3 we give a concrete formula 
describing these weight multiplicities as the number of integer points in certain rational convex 
polytopes; we indicate that this is again amenable to algorithmic implementation. Finally, we 
show that in the limit we recover the corresponding statements of |§4.2| 

6.1. The Semi-Classical Limit. Since the JC-action on M C P{V) originates from a linear 
action on V, each graded part of the homogeneous coordinate ring C[M] is naturally a finite- 
dimensional ii'-representation and can thus be decomposed into irreducible sub-representations, 

oo 

C[M]^0 Vx0ilomK{VxX[M]k). 
k=o AeA*nt:^ 

We shall encode their multiplicities, suitably re-scaled, in the following sequence of discrete 
measures, 

(6.1) f^M,k--=^ E dimHom;,(y;,C[M],)5,/fc. 

AeA*nt^ 

The factor accommodates for the growth of the dimension of a generic irreducible representation 
in the coordinate ring, which has highest weight in tt,Q (we still assume that Assumption 1.1 is 
in place). 

It is well-known that in the semi-classical limit A: — )• oo this sequence of measures converges in 
distribution to the non-Abelian Duistermaat-Heckman measure [501 SSl iHSl [79l IHO l 1102] . 

(6.2) /.iS^.^DHf,. 

In fact, one can show using the Hirzebruch-Riemann-Roch theorem that the piecewise polynomial 
density function of the Duistermaat-Heckman measure is at any regular point A G tt.g for 
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(equivalently, for ^t) given by 

(6.3) lim ^ dimRomK{Vkx*,C[M]k). 

Here, D is the number of elements in the generic stabilizer of points in Mj^and the right-hand 
side is understood as the limit over the subsequence with k\ integral. The additional factor k"^ 
comes from the fact that we now consider the density with respect to the r-dimensional Lebesgue 
measure d\. 



It is well-known that the right-hand side multiplicity in (6.3) is a quasi-polynomial in k, i.e., 
a polynomial whose coefficients are periodic functions of k. Observe that the existence of the 
limit implies that (a) the degree of this quasi-polynomial is at most n — R — r, and (b) if it 
is of maximal degree then its leading order coefficient is in fact a constant independent of k 
(namely, the limit). We shall therefore call the limit the maximal- order growth coefficient of the 
quasi-polynomial. 

The rational points of the moment polytope Ak{M) are precisely those of the form X/k with 
C C[M]k IMJ. In other words, suppDH^nt^ = UfoSupp/i^ ^. 

Restricting the action to the maximal torus T C if, we observe that the Abelian Duistermaat- 
Heckman measure DH^ captures the asymptotic distribution of weights in the homogeneous 
coordinate ring of M, i.e., the asymptotics of the character of C[Af]fc as A; — )• oo. 

Example 6.1. For strictly dominant and integral A G A* n tJ^g, the Borel-Weil theorem shows 
that the homogeneous coordinate ring of the coadjoint orbit Ox C PiVx) is equal to 

oo 

C[Ox] = ^V,\. 

k=0 

Therefore, all the multiplicity measures /i^^ ^ (and hence their limit) are equal to the Dirac 
measure at A. 

6.2. Multiplicities of Irreducible Representations via Finite Differences. Multiplicities 
of weights and highest weights in finite-dimensional if-representations are related by iteratively 
taking (negative) finite differences in the directions of the positive roots. This can be seen as the 



"quantized" version of the derivative principle (Corollary 3.4). Its proof is in essence a rephrasing 



of the Weyl character formula, an idea which goes back at least to Steinberg |97j . 

Lemma 6.2. Denote by nix and tjit the highest weight and weight multiplicity function, re- 
spectively, of a finite-dimensional K -representation V. Then on the positive Weyl chamber we 
have 



rriK = JJ -Da rriT 



\a>0 / 

■+ 

where (Dam)(X) = m(A -\- a) — m(A) is the finite- difference operator in direction a. Note that 
any two of the operators Da commute, so that their product is independent of the order of 
multiplication. 

Proof By linearity of the finite-difference operators it suffices to establish the lemma for a single 
irreducible representation Vx of highest weight A. It will be convenient to work with the formal 
character ch(yx) = mj'(/i) |19| I65| . By the Weyl character formula, 

]J (1 - e"") ch(yx) = e-P (-l)'("')e"'(^+^) 

where W is the Weyl group, l{w) the length of a Weyl group element w, and p half the sum of 
the positive roots. 

Now observe that the left-hand side is the generating function of (nQ>o ~Da) "^T, since 
taking finite differences corresponds to multiplying the generating function by 1 — e~". Up to 



Quotienting out a discrete subgroup leaves the Duistermaat-Heckman measure invariant, but changes the 
weight lattice, and therefore the normalization of the Lebesgue measure d\. 
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terms corresponding to non-dominant weights, the right-hand side is equal to e^, 
generating function of rriK = ^x- ■ The assertion follows from this. 



which is the 
□ 



Note that the Weyl character formula can be seen as the representation-theoretic analogue 
of the Harish-Chandra formula that was used to establish [Corollary 3.4 In the semi-classical 
limit (6.2), the finite differences become infinitesimal and we recover an alternative proof of 



Corollary 3.4 in the algebro-geometric setting. 



This argument can also be turned aro und to establis h (6.2) for general compact Lie groups K 
from its Abelian version [44, (34.8)] and [Corollary 3.4|p " 



6.3. Multiplicities for Projective Spaces. As in §4.2 let M = F{V) be the complex projec- 
tive space for a unitary /C-representation V. Its homogeneous coordinate ring is equal to the 
symmetric algebra, 

oo 

C[¥{V)] = Sym{V) = ^Sym\V). 

k=0 

Choose a weight-space decomposition V = ^^^qCv^, and identify V = C""^^ and U(y) = 
U(n -|- 1) accordingly. Observe that the maximal torus T K acts via the standard maximal 
torus of U(n -|- 1), that is, the set of unitary diagonal matrices, which we denote by T. 

Each symmetric tensor power Sym'^(C""*'^) is an irreducible representation of U(n -|- 1). 
Its weight spaces are all one- dimensional, and the weights that occur are precisely the A = 
diag(Ao, • • • , A„) with i\j E Z, \j > 0, and \j = k |33| . Clearly, we can identify this set of 
weights with the integral points in A;A„, where A„ is the n-dimensional standard simplex in M"^-'^. 
In the language of Young diagrams, these are the weight vectors corresponding to semistandard 
Young tableaux of shape {k) with entries in {0, . . . , n}. 

To determine the weight multiplicities with respect to T C i^, we have to "restrict" each 
weight to t. This corresponds precisely to applying the map P: W"^"^ — )• t*,(tfc) i— )• X^^tfcWfc 
introduced in Proposition 4.10 Therefore, the multiplicity in Sym^iV) of a weight A G A* is 



given by counting integral points in a rational convex polytope parametrized by k and A: 

(6.4) mT,fc(A) = #(A(A,A:)nZ"+i), 

where 



(6.5) A(A, k) = |(tj) G M"+^ : tj > 0, ^ tjUjj = \^tj = k\ 

j=0 3=0 

Remark 6.3. Such vector partition functions can be evaluated efficiently using Barvinok's algorithm 
if the group K and the ambient dimension dim y = n -|- 1 is fixed (H El [6] , namely in time 
0(poly(log /;:)). In fact, rriT^k is a piecewise quasi-polynomial function in both A and k, and there 
are parametric generalizations of Barvinok's algorithm for computing these quasi-polynomials 
[ion llOOj . Since we can compute multiplicities of irreducible ^-representations by taking finite 



differences of weight multiplicities in the direction of positive roots (Lemma 6.2), this can also be 



4.2 can be used to 



done efficiently if K is fixed. We will report on a generalization of this technique to the general 
branching problem for compact connected Lie groups in a forthcoming article [21j. 

There is also a jump formula by Boysal and Vergne [12j, which as in 
inductively compute the quasi-polynomials chamber by chamber. 

We now turn to the semi-classical limit. As A; — t- oo, it is clear that 

1 



AeA„niz"+i 



converges to Lebesgue measure on the standard simplex A„, normalized to total volume 



hm - — dim Sym [V ] = iim - — 



1 



'^We thank Allen Knutson for pointing tliis out, as well as for sketching a self-contained proof of (|6.2[) 
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Therefore, /^p^^y-j converges to the push-forward of Lebesgue measure on A„ along the map P. 



By the semi-classical limit (6.2), this is of course equivalent to the assertion of Proposition 4.10 
Moreover, note that the quantity 

mT,k{k\) = # (A(fcA, k) n = # (A(A, 1) n 1^"+^) 

is the Ehrhart quasi-polynomial associated to rational polytope A(A, 1) [8j. It is intuitively clear 
that its growth in k should be related to the volume of this polytope. Indeed, 

mT,k{k\) = A;"-WolA(A,l) + 0(A;"-''-i), 



where vol is the (n — r)-dimensional volume with respect to the measure dt/dX defined in 



Exercise 3.29]. Observe that this agrees with (6.3) and Proposition 4.12} The maximal-order 



(4.2 



growth coefficient is a constant equal to the Abelian Duistermaat-Heckman density at point A. 

7. KRONECKER AND PLETHYSM COEFFICIENTS 

In this section, we describe the representation theory of the quantum marginal problem in 
more detail. For distinguishable particles, the relevant multiplicities can be expressed in terms of 



decomposing tensor products of irreducible representations of the symmetric group (§7.1 ). In 



particular, the joint eigenvalue distribution of the reduced density matrices of a tripartite pure 



state is determined by the asymptotics of the Kronecker coefficients (see (7.4 )). We emphasize that 



by specializing the method described in §6.3 we get a novel algorithm for computing Kronecker 



coefficients which is efficient for Young diagrams of bounded height. Indistinguishable particles 



correspond to certain plethysm coefficients and we conclude by illustrating this connection (§7.2). 



7.1. Kronecker Coefficients. Recall that for N distinguishable particles we have to consider 
the action of K = S\]{di) x . . . x SU((ijv) on a coadjoint SU((ii • • • dAf)-orbit M = O^, where we 
now assume that A is an integral weight in tI,Q. The multiplicity measures /i^^ are determined 
by the decomposition of the homogeneous coordinate ring 



d\---dj^ 

kX 
k=0 



into JC-isotypical components (the superscript labels the corresponding SU). 

We can express this equivalently using the representation theory of the symmetric group Sm- 
Recall that by Schur-Weyl duality the diagonal action of SU(d) and the permutation action of 
Sm on (C^)®™ generate each other's commutant, so that 

(7.1) (C'^)«™^0y^rf^[^]. 

Here, the sum runs over all Young diagrams fi = (/xi, . . . , /x^) with |/i| := Hj = m boxes and at 
most d rows, is the irreducible representation of SU(d) with highest weight X i-^ i Xjfij, 
and [^] is the corresponding irreducible representation of Sm (see |33| for details). We shall freely 
identify Young diagrams and the corresponding highest weights. 

In particular, we can realize the irreducible representation in (C*^^ . . . (S> C'^N^^\kX\ ^ 



Comparing the Schur-Weyl decomposition (7.1 ) for the full Hilbert space with the tensor product 



of the decompositions for the individual subsystems, we find that 

V^y-'^= <®...®y,'^;®Hom5,^,|([fcA],[Ai]0...0[A^]), 

Ai,...,Ajv 

where the sum runs over the Young diagrams Aj with |fcA| boxes and at most di rows. Therefore, 
/^M,fc =]J^^ dimHomS|^-|([A:A], [Ai]* . . . [Aa^]*) 5(Ai/fc,...,Ajv/fc) 

(7.2) 



1 



Ai,...,Ajv 



dimHoms ([/cA], [Ai] ® [Ajv]) (^(Ai/fc,...,Ajv/fc), 
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Figure 7. (a) and (b) Illustration of the multiplicity measures for the 

mixed-state quantum marginal problem of two qubits with global spectrum 
(4/7,2/7,1/7,0) and k = 28,56, which have been computed by the algorithm 
described in §7.1| (c) Their semi-classical limit, i.e., the corresponding Duister- 
maat-Heckman measure as computed in §5.2 



where the latter identity holds due to the self-duality of the representations of the symmetric 
group. In particular, the rational points of the non-Abelian moment polytope /S.K{Oj^ are 
precisely 



(7.3) 



IJ (Ai/A:,...,A^/A:): 



[k\] C [Ai] (g) ... (g) [A at], where the 

Aj have |A;A| boxes and at most di rows 



See Figure 7 for an illustration of the multiplicity measures corresponding to the mixed-state 



quantum marginal problem for two qubits discussed in §5.2 



Remark 7.1. We can write the multiplicities in (7.2) in the following symmetric form: 

■ ^' 

dimHoms, -,([A;A]> [^i]* [Xn]*) = dim [Ai]* (g) . . . (g [Xn]* ® [kX]* 



\k\\ 



Observe that the right-hand side is a multiplicity for the pure-state quantum marginal problem 
for (g) . . . (g C'^'^ (g C'^^-'^N _ Indeed, the homo geneous coordinate ring of a projective space is 



just the symmetric algebra (§6.3), whose graded parts correspond to the trivial representations of 



the symmetric groups. This is the representation-theoretic perspective on purification (cf. §2.3 
in particular Proposition 2.8). 



For the tripartite pure-state quantum marginal problem (equivalently, the mixed-state bipartite 
quantum marginal problem), the relevant multiplicities are the well-known Kronecker coefficients 
of the symmetric group, 

gx,^,,u = dim{[X]0[fi](g)[iy]fK 

They are the symmetric group analogue of the Littlewood-Richardson coefficients of the unitary 
group (in fact, the latter can be considered as a special case) but much harder to compute in 
general, since there is no combinatorial description like the Littlewood-Richardson rule. 



The corresponding characterization (7.3) of the non-Abelian moment polytope has already 



been observed in |23| 1611 [22], as well as in j28| for the projection onto two of the subsystems. 
The semi-classical limit refines this characterization: Not only can one read off the existence of 
quantum states with given marginal eigenvalue spectra from the asymptotic non- vanishing of 
the corresponding Kronecker coefficients gk\,kfi,ku, but their growth also encodes the probability 
of finding these eigenvalue spectra when the global state is chosen according to the invariant 
probability measure. Explicitly, ( |6.2[ ) states that 

1 



(7.4) 



kP ^ 



9X,fi,iy ^\/k,n/k,v/k DH 





SU{a)xSU(6)xSU(c)' 



where p = n — R = abc — 1 — a{a — l)/2 -|- b{b — l)/2 -|- c(c — l)/2, and where the sum runs over 
all Young diagrams A, /x, with k boxes and at most a, b and c rows, respectively. 



The method described in §6.3 in particular provides a novel algorithm for computing the 



Kronecker coefficients which is efficient for Young diagrams of bounded height: Using the finite- 



difference formula of Lemma 6.2 we can reduce to the computation of a bounded number of 
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weight multiplicities (6.4), which using Barvinok's algorithm can be evaluated in polynomial 



time in the input size, i.e., in time 0(poly(log fc)), where k is the number of boxes of the Young 



diagrams. As mentioned in Remark 6.3[ we will elaborate on this algorithm in a forthcoming 
article l211. 



7.2. Plethysm Coefficients. While the quantum marginal problem for distinguishable particles 
is connected to (generalized) Kronecker coefhcients, it is for indistinguishable particles related to 
certain plethysm coefficients. Indeed, if M = P(Va) for an irreducible SU((i)-representation V\ 
then its coordinate ring consists of the plethysms 

oo 

C[M] = 0Sym^(yA). 

A;=0 

See e.g. |78| for more information on plethysms, which are in general defined as the composition 
of Schur functors. In particular, the bosonic and fermionic pure-state marginal problem are 
related to the asymptotics of Sym'^(Sym^(C'^)) and Sym'^(A^(C'^)), respectively, as /c — ?■ oo. 

Let us illustrate this by giving an alternative derivation of the Duistermaat-Heckman measures 



for N bosonic qubits (cf. §5.3). We will explicitly compute the asymptotic weight multiplicity 
distribution of the plethysm Sym'^(Sym'^(C^)) as /c — )• oo, and then apply the derivative principle. 
The main combinatorial tool we shall employ are the g-binomial coefficients 



n 



Recall that these are defined in terms of the g-integers [n\q 



n „ n ■ 



and (7- factorials [n]q\ 

N I 



l]q . . . [V\q. We start with the following description of the character of Sym (Sym (C )) 



in terms of g- binomial coefficients: 



ch ^Sym^ (Sym 



N 



and q - 


= e'^i . : 




'k + N' 


)))- 


N 



-.-kN 



Proposition 7.3. As functions on the open unit disk {</ € C : |g| < 1} one has for fixed iV G N 
and A; —7- 00 the following asymptotic equivalence 

'k + N' 

N 



[k] 



Proof. Following [55, (9.1)], for any fixed c G N and \q\ < 1 one has 

1 - 1 
lim = . 

fc^oo 1 — q ^ — q 

By applying this identity both to the numerator and the denominator, 

[k + C]q 



lim 

fc— ^00 



1. 



Hence 



lim 

A;— >oo 



lim 



[k + N]q[k + N - l]q . 



\k + l] 



□ 



The following corollary is an easy application of Osgood's theorem f85] to Proposition 7.3| (see 
e.g. [71)- 



Corollary 7.4. Fix N, and define fk{q) = ['^a/^] g ^^'^ 9k{q) = [k]^ ■ Then, the sequences 
fk and gi^ converge, as k ^ 00, pointwise to the same holomorphic function on some open dense 
subset in the open unit disk {q & C : \q\ < 1}. In particular, the limits are equal as power series. 



We can use this result to extract asymptotic multiplicity information. 
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Proposition 7.5. The discrete measures := /^Sym^(c2) ^, ^■^ defined in (6.1) tend in the 



U{1) 



limit k 



oo to 



times the probability distribution of the sum of N independent random 



variables uniformly distributed on 

U(l) 



Proof. By 



Proposition 7.2 



{N),k 



[-1,1]- 

is a finite measure with generating function 



X J U(l) / N 



1 



'k + N' 
N 



Let Vk, ojk be finite measures with generating functions [A^]„2/fc! and 



1 

w 



N 



k 




respectively. Obviously, uj^ is asymptotically distributed like the sum of inde pendent random 

U(l) 



variables uniformly distributed on the interval [—1,1], and by 
Since Uk N\ 6q as k ^ oo, this implies our assertion. 



Corollary 7.4 



so IS Ai(7^);fc * ^k- 

□ 



By using the semi-classical limit (6.2), we conclude once again that the Abelian Duister- 
maat-Heckman measure is given by the formula that was established in Proposition 5.11 The 
non- Abelian Duistermaat-Heckman measure is obtained as in Corollary 5.12 by applying the 
derivative principle. 

Example 7.6 (A^ = 2). For the plethysms Sym'^(Sym^(C^)) we can also illustrate the semi-classical 
limit for the SU(2)-action, since the decomposition into irreducible SU(2)-representations is 
well-known |78| §1.5, Example 6 (a)]: 

Sym*^(Sym2(C2)) ^ Sym2^(C2) Sym^^-^CC^) ... Sym^'^iC'^) 

The last summand is Sym^(C^) for odd /c, and Sym^(C^) for even k. Therefore, the discrete 
measures as defined in (6.1) are given by 

1 . 



SU(2) 

P(Sym2(C2)),fc 



=2fc,2fc-4,. 



,2|0 



A;2- 



In the limit /c — )■ oo, they converge to the non- Abelian Duistermaat-Heckman measure as 
computed in Corollary 5.12 



dh: 



SU(2) 

P(Sym2(C2)) 



■\ {{x + 2)0 - 2x1 + 2)+) dx = |l[o,2)(x)dx. 



See Figure 6 for an illustration. 

Remark 7.7. The description of the character of Sym^(Sym^(C^)) via g-binomial coefficients has 
the additional advantage that one is able to compute all higher cumulants and moments of the 
associated distribution for any fixed k and A^ (see [90]). This is due to a method by Panny 



Acknowledgements 

We would like to thank Alonso Botero, Emmanuel Briand, Peter Biirgisser, Benoit Collins, 
David Gross, Christian Ikenmeyer, Markus P. Miiller, Mercedes Rosas, and Volkher Scholz for 
helpful discussions. We thank Graeme Mitchison for joint initial discussions on the topic of 
asymptotics of Kronecker coefficients. The second author would like to express his particular 
gratitude to Frances Kirwan for many fruitful discussions regarding moment maps and invariant 
theory and their many uses. 

This work is supported by the Swiss National Science Foundation (grants PP00P2-128455 and 
200021_138071), the German Science Foundation (SFB/TR12, and grants CH 843/1-1 and CH 
843/2-1), the National Center of Competence in Research 'Quantum Science and Technology', 
and the Excellence Initiative of the German Federal and State Governments through the Junior 
Research Group Program within the Institutional Strategy ZUK 43. 



40 m. christandl, b. doran, s. kousidis, and m. walter 

References 

1. G. Aubmn, S. Szarek, and E. Werner, HasUng's Additivity Counterexample via Dvoretzky's Theorem, 
Communications in Mathematical Physics 305 (2011), 85-97. 

2. G. Aubrun, S. Szarek, and D. Ye, Entanglement thresholds for random induced states, arXiv: 1106.2264, 2011. 

3. , Phase transitions for random states and a semi-circle law for the partial transpose, arXiv: 11 12.4582, 

2011. 

4. A. Barvinok, Computing the Volume, Counting Integral Points, and Exponential Sums, Discrete and 
Gomputational Geometry 10 (1993), 123-141. 

5. , A Polynomial Time Algorithm for Counting Integral Points in Polyhedra when the Dimension is 

Fixed, Mathematics of Operations Research 19 (1994), 769-779. 

6. A. Barvinok and J. E. Pommersheim, New Perspectives in Algebraic Combinatorics, MSRI Publications, 
vol. 38, ch. An Algorithmic Theory of Lattice Points in Polyhedra, Cambridge University Press, 1999. 

7. A. F. Beardon and D. Minda, On the pointwise limit of complex analytic functions. The American Mathe- 
matical Monthly 110 (2003), 289-297. 

8. M. Beck and S. Robins, Computing the Continuous Discretely: Integer-Point Enumeration in Polyhedra, 
Springer, 2009. 

9. A. Berenstein and R. Sjamaar, Coadjoint orbits, moment polytopes, and the Hilbert-Mumford criterion. 
Journal of the American Mathematical Society 13 (2000), 433-466. 

10. N. Berline, E. Getzler, and M. Vergne, Heat Kernels and Dirac Operators, Springer- Verlag, 2003. 

11. R. E. Borland and K. Dennis, The conditions on the one-matrix for three-body fermion wavefunctions with 
one-rank equal to six, Journal of Physics B 5 (1972), 7-15. 

12. A. Boysal and M. Vergne, Paradan's wall crossing formula for partition functions and Khovanski-Pukhlikov 
differential operators, Annales de ITnstitut Fourier 59 (2009), 1715-1752. 

13. S. Bravyi, Requirements for compatibility between local and multipartite quantum states. Quantum Information 
and Computation 4 (2004), 012-026. 

14. M. Brion, Sur I'image de I'application moment, Seminaire d'Algebre Paul Dubreil et Marie-Paule Malliavin, 
Lecture Notes in Mathematics, vol. 1296, Springer, 1987, pp. 177-192. 

15. P. Biirgisser, M. Christandl, and C. Ikenmeyer, Even Partitions in Plethysms, Journal of Algebra 328 (2011), 
322-329. 

16. , Nonvanishing of Kronecker coefficients for rectangular shapes. Advances in Mathematics 227 (2011), 

2082-2091. 

17. P. Biirgisser, J. M. Landsberg, L. Manivel, and J. Weyman, An Overview of Mathematical Issues arising 
in the Geometric Complexity Theory approach to VP =^ VNP, SIAM Journal on Computing 40 (2011), 

1179-1209. 

18. A. Cannas da Silva, Lectures on Symplectic Geometry, corr., 2nd printing ed.. Lecture Notes in Mathematics, 
vol. 1764, Springer- Verlag, 2008. 

19. R. W. Carter, G. Segal, and I. G. MacDonald, Lectures on Lie groups and Lie algebras, London Mathematical 
Society, 1995. 

20. M. Christandl, A Quantum Information- Theoretic Proof of the Relation between Horn's Problem and the 
Littlewood-Richardson Coefficients, Proceedings of Computability in Europe: CiE 2008, Lecture Notes in 
Computer Science, vol. 5028, Springer, 2008, pp. 120-128. 

21. M. Christandl, B. Doran, and M. Walter, Computing Multiplicities of Lie Group Representations, in 
preparation, 2012. 

22. M. Christandl, A. W. Harrow, and G. Mitchison, On Nonzero Kronecker Coefficients and their Consequences 
for Spectra, Communications in Mathematical Physics 270 (2007), 575-585. 

23. M. Christandl and G. Mitchison, The Spectra of Quantum States and the Kronecker Coefficients of the 
Symmetric Group, Communications in Mathematical Physics 261 (2006), 789-797. 

24. M. Christandl and A. Winter, "Squashed Entanglement" - An Additive Entanglement Measure, Journal of 
Mathematical Physics 45 (2004), 829-840. 

25. A. J. Coleman, Structure of Fermion Density Matrices, Reviews in Modern Physics 35 (1963), 668-686. 

26. A. J. Coleman and V. I. Yukalov, Reduced Density Matrices: Coulson's Challenge, Lecture Notes in Chemistry, 
vol. 72, Springer, 2000. 

27. B. Collins, I. Necliita, and D. Ye, The absolute positive partial transpose property for random induced states, 
arXiv:1108.1935, 2011. 

28. S. Daftuar and P. Hayden, Quantum state transformations and the Schubert calculus. Annals of Physics 315 

(2004), 80-122. 

29. A. H. Dooley, J. Repka, and N. J. Wildberger, Sums of adjoint orbits. Linear and Multilinear Algebra 36 
(1993), 79-101. 

30. J. J. Duistermaat and G. J. Heckman, On the Variation in the Cohomology of the Symplectic Form of the 
Reduced Phase Space, Inventiones Mathematicae 69 (1982), 259-268. 

31. W. Feller, An introduction to probability theory and its applications. Vol. II., Second edition, John Wiley & 
Sons Inc., 1971. 



EIGENVALUE DISTRIBUTIONS OF REDUCED DENSITY MATRICES 



41 



32. A. Frumkin and A. Goldbcrgcr, On the distribution of the spectrum of the sum of two hermitian or real 
symmetric matrices, Advances in Applied Mathematics 37 (2006), 268-286. 

33. W. Pulton, Young Tableaux, Student Texts, London Mathematical Society, 1997. 

34. , Eigenvalues, invariant factors, highest weights, and Schubert calculus, Bulletin of the American 

Mathematical Society 37 (2000), 209-249. 

35. S. Goldstein, J. L. Lebowitz, R. Tumulka, and N. Zanhi, Canonical typicality, Physical Review Letters 96 
(2006), 050403. 

36. D. M. Greenberger, M. A. Horne, and A. Zeilinger, Going beyond Bell's Theorem, Bell's Theorem, Quantum 

Theory, and Conceptions of the Universe (M. Kafatos, ed.), Khiwor. 1089, pp. 69-72. 

37. V. Guillemin, E. Lcrman, and S. Sternberg, On the Kostant multiplicity formula, Journal of Geometry and 
Physics 5 (1988), 721-750. 

38. , Symplectic Fibrations and Multiplicity Diagrams, Cambridge University Press, 1996. 

39. V. Guillemin and E. Prato, Heckman, Kostant, and Steinberg formulas for symplectic manifolds, Advances in 
Mathematics 82 (1990), 160-179. 

40. V. Guillemin and R. Sjamaar, Convexity Properties of Hamiltonian Group Actions, American Mathematical 
Society, 2005. 

41. V. Guillemin and S. Sternberg, Geometric Asymptotics, revised edition ed.. Mathematical Surveys and 
Monographs, vol. 14, American Mathematical Society, 1977. 

42. , Convexity Properties of the Moment Mapping, Inventioucs Matlicmaticac 67 (1982), 491 513. 

43. , Geometric Quantization and Multiplicities of Group Representations, Inventiones Mathematicae 67 

(1982), 515-538. 

44. , Symplectic Techniques in Physics, Cambridge University Press, 1984. 

45. Harish-Chandra, Differential operators on a semisimple Lie algebra, American Journal of Mathematics 79 
(1957), 87-120. 

46. M. B. Hastings, Superadditivity of communication capacity using entangled inputs, Nature Physics 5 (2009), 
255-257. 

47. P. Hayden, D. Leung, P. W. Shor, and A. Winter, Randomizing Quantum States: Constructions and 
Applications, Communications in Mathematical Physics 250 (2004), 371-391. 

48. P. Hayden, D. W. Leung, and A. Winter, Aspects of generic entanglement. Communications in Mathematical 
Physics 265 (2006), 95-117. 

49. P. Hayden and J. Preskill, Black holes as mirrors: quantum information in random subsystems, Journal of 
High Energy Physics 2007 (2007), 120. 

50. G. J. Heckman, Projections of Orbits and Asymptotic Behaviour of Multiplicities for Compact Connected Lie 
Groups, Inventiones Mathematicae 67 (1982), 333-356. 

51. U. Helmke and J. Rosenthal, Eigenvalue inequalities and Schubert calculus, Mathematische Nachrichten 171 
(1995), 207-225. 

52. A. Higuchi, On the one-particle reduced density matrix of a pure three-qutrit quantum state, arXiv:quant- 

ph/0309186v2, 2003. 

53. A. Higuchi, A. Sudbery, and J. Szulc, One-qubit reduced states of a pure many-qubit state: polygon inequalities, 
Physical Review Letters 90 (2003), 107902. 

54. K. Huang, Statistical Mechanics, Wiley, John & Sons, 1990. 

55. V. Kac and P. Cheung, Quantum Calculus, Univcrsitcxt, Springer- Vcrlag, 2002. 

56. A. A. Kirillov, Merits and demerits of the orbit method, Bulletin of the American Mathematical Society 36 
(1999), 433-488. 

57. A. Kirillov Jr., An Introduction to Lie Groups and Lie Algebras, Cambridge studies in advanced mathematics, 
Cambridge University Press, 2008. 

58. F. Kirwan, Cohomology of quotients in symplectic and algebraic geometry, Mathematical Notes, Princeton 
University Press, 1984. 

59. , Convexity properties of the moment mapping. III, Inventiones Mathematicae 77 (1984), 547-552. 

60. A. Klyachko, Stable vector bundles and Hermitian operators, Selecta Mathematica, New Series 4 (1998), 
419-445. 

61. , Quantum marginal problem and representations of the symmetric group, arXiv:quant-ph/0409113, 

2004. 

62. , Quantum marginal problem and N-representability, Journal of Physics: Conference Series 36 (2006), 

72-86. 

63. , The Pauli exclusion principle and beyond, arXiv:0904.2009, 2009. 

64. A. Klyachko and M. Altunbulak, The Pauli Principle Revisited, Communications in Mathematical Physics 
282 (2008), 287-322. 

65. A. Knapp, Lie Groups: Beyond an Introduction, second edition ed.. Progress in Mathematics, vol. 140, 

Birkhauscr, 2002. 

66. A. Knutson, The symplectic and algebraic geometry of Horn's problem, Linear Algebra and its Applications 
319 (2000), 61-81. 



42 



M. CHRISTANDL, B. DORAN, S. KOUSIDIS, AND M. WALTER 



67. , Schuhert Calculus and Quantum Information, Quantum Marginals and Density Matrices Workshop, 

Field Institute, Toronto, 2009. 

68. A. Knutson and T. Tao, The honeycomb model o/GL„(C) tensor products I: Proof of the saturation conjecture, 
Journal of the American Mathematical Society 12 (1999), 1055-1090. 

69. , Honeycombs and Sums of Hermitian Matrices, Notices of the American Mathematical Society 38 

(2001), 175-186. 

70. A. Knutson, T. Tao, and C. Woodward, The honeycomb model of GL„{C) tensor products II: Puzzles 
determine facets of the Littlewood-Richardson cone. Journal of the American Mathematical Society 17 (2003), 
19-48. 

71. E. Lerman, E. Meinrcnkcn, S. Tolman, and C. Woodward, Non-abehan convexity by symplectic cuts. Topology 
37 (1998), 245-259. 

72. B. V. Lidskii, Spectral polyhedron of a sum of two Hermitian matrices, Functional Analysis and Its Applications 
16 (1982), 139-140. 

73. Y.-K. Liu, Consistency of Local Density Matrices is QMA-complete, Proceedings of RANDOM 2006, pp. 438- 
449. 

74. Y.-K. Liu, M. Christandl, and F. Verstraete, Quantum Computational Complexity of the N -Representability 
Problem: QMA Complete, Physical Review Letters 98 (2007), 110503. 

75. S. Lloyd, Excuse our ignorance. Nature Physics 2 (2006), 727-728. 

76. S. Lloyd and H. Pagels, Complexity as thermodynamic depth. Annals of Physics 188 (1988), 186-213. 

77. E. Lubkin, Entropy of an n-system from its correlation with a k-reservoir, Journal of Mathematical Physics 
19 (1978), 1028-1031. 

78. I. G. MacDonald, Symmetric Functions and Hall Polynomials, Oxford Mathematical Monographs, 1995. 

79. E. Meinrenken, On Riemann-Roch Formulas for Multiplicities, Journal of the American Mathematical 
Society 9 (1996), 373-389. 

80. E. Meinrenken and R. Sjamaar, Singular reduction and quantization. Topology 38 (1999), 699-762. 

81. M. P. Miiller, O. C. O. Dahlsten, and V. Vedral, Unifying typical entanglement and coin tossing: on 
randomization in probabilistic theories, arXiv:1107.6029, 2011. 

82. K. Mulmuley, Geometric complexity theory VI: the flip via saturated and positive integer programming 
in representation theory and algebraic geometry. Technical report. Computer Science Department, The 
University of Chicago, 2007. 

83. K. Mulmuley and M. Sohoni, Geometric complexity theory I: An approach to the P vs. NP and related 
problems, SIAM Journal on Computing 31 (2001), 496-526. 

84. , Geometric complexity theory IT. Towards explicit obstructions for embeddings among class varieties, 

SIAM Journal on Computing 38 (2008), 1175-1206. 

85. W. F. Osgood, Note on the functions defined by infinite series whose terms are analytic functions of a 
complex variable; with corresponding theorems for definite integrals, Annals of Mathematics. Second Series 3 
(1901/02), 25-34. 

86. D. N. Pago, Average entropy of a subsystem. Physical Review Letters 71 (1993), 1291-1294. 

87. , Black hole information. Proceedings of the 5th Canadian Conference on General Relativity and 

Relativistic Astrophysics, American Mathematical Society, 1994. 

88. W. Panny, A note on the higher moments of the expected behavior of straight insertion sort, Information 
Processing Letters 22 (1986), 175 177. 

89. S. Popescu, A. J. Short, and A. Winter, Entanglement and the foundations of statistical mechanics. Nature 
Physics 2 (2006), 754-758. 

90. H. Prodinger, On the moments of a distribution defined by the Gaussian polynomials. Journal of Statistical 
Planning and Inference 119 (2004), 237-239. 

91. N. Rcssayre, Geometric invariant theory and the generalized eigenvalue problem, Inventiones Mathematicae 
180 (2010), 389-441. 

92. M. B. Ruskai, N -Representability Problem: Conditions on Geminals, Physical Review 183 (1969), 129-141. 

93. , Connecting N -representability to Weyl's problem: the one-particle density matrix forn = 3 and 

r = 6, Journal of Physics A 40 (2007), F961-F967. 

94. P. W. Shor, Equivalence of Additivity Questions in Quantum Information Theory, Communications in 
Mathematical Physics 246 (2004), 453-472. 

95. R. Sjamaar, Holomorphic slices, symplectic reduction and multiplicities of representations. Annals of Mathe- 
matics 141 (1995), 87-129. 

96. T. A. Springer, Invariant theory. Lecture Notes in Mathematics, Vol. 585, Springer- Verlag, 1977. 

97. R. Steinberg, A general Clebsch-Gordan theorem. Bulletin of the American Mathematical Society 67 (1961), 
406-407. 

98. F. H. Stillinger, Mathematical challenges from theoretical/ computational chemistry. National Academy Press, 
1995. 

99. L. C. Vcnuti and P. Zanardi, Probability density of quantum expectation values, arXiv: 1202.4810, 2012. 
100. S. Verdoolaogc and M. Bruynooghe, Algorithms for Weighted Counting over Parametric Polytopes: A Survey 

and a Practical Comparison, ITSL 2008 (2008), 60-66. 



EIGENVALUE DISTRIBUTIONS OF REDUCED DENSITY MATRICES 



43 



101. S. Verdoolaege, R. Seghir, K. Beyls, V. Loechner, and M. Bruynooghe, Counting Integer Points in Parametric 
Polytopes using Barvmok's Rational Functions, Algorithmica 48 (2007), 37-66. 

102. M. Vergne, Quantization of Algebraic Cones and Vogan's Conjecture, Pacific Journal of Mathematics 182 
(1998), 113-135. 

103. N. M. J. Woodliouse, Geometric quantization, second edition ed., Oxford Mathematical Monographs, The 
Clarendon Press, Oxford University Press, 1992. 

104. C. T. Woodward, Localization for the norm-square of the moment map and the two-dimensional Yang-Mills 
integral. Journal of Symplectic Geometry 3 (2005), 17-54. 

105. K. Zyczkowski and H.-J. Sommers, Induced measures in the space of mixed quantum states, Journal of Physics 
A 34 (2001), 7111-7125. 

Matthias Christandl, Institute for Theoretical Physics, ETH Zurich, Wolfgang-Pauli-Strasse 
27, CH-8093 Zurich, Switzerland 

E-mail address: |christand l@phys . ethz . ch| 

Brent Doran, Department of Mathematics, ETH Zurich, Ramistrasse 101, CH-8092 Zurich, 
Switzerland 

E-mail address: Ibrent . doranSmath . ethz . chl 

Stavros Kousidis, Institute for Theoretical Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, 
CH-8093 Zurich, Switzerland 

Institute of Physics, University of Freiburg, Rheinstrasse 10, 79104 Freiburg, Germany 
E-mail address: |st . kousidis@googlemail . com| 

Michael Walter, Institute for Theoretical Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, 
CH-8093 Zurich, Switzerland 

E-mail address: mwalt erOitp . phys . ethz . ch| 



